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Abstract 

The spectacular experimental results of the last few years in cavity 
quantum electrodynamics and trapped ions research has led to very high 
level laboratory performances. Such a stimulating situation essentially 
stems from two decisive advancements. The first is the invention of re- 
liable protocols for the manipulation of single atoms. The second is the 
ability to produce desired bosonic environments on demand. These pro- 
gresses have led to the possibility of controlling the form of the coupling 
between individual atoms and an arbitrary number of bosonic modes. As 
a consequence, fundamental matter-radiation interaction models like, for 
instance, the JC model and most of its numerous nonlinear multiphoton 
generalizations, have been realized or simulated in laboratory and their 
dynamical features have been tested more or less in detail. This topical 
paper reviews the state of the art of the theoretical investigations and of 
the experimental observations concerning the dynamical features of the 
coupling between single few-level atoms and two bosonic modes. In the 
course of the paper we show that such a configuration provides an excellent 
platform for investigating various quantum intermode correlation effects 
tested or testable in the cavity quantum electrodynamics and trapped ion 
experimental realms. In particular we discuss a mode-mode correlation 
effect appearing in the dynamics of a two-level atom quadratically cou- 
pled to two bosonic modes. This effect, named parity effect, consists in a 
high sensitivity to the evenness or oddness of the total number of bosonic 
excitations. 
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1 Foreword 



The history of the basic ideas of quantum theory, is studded with proposal of 
gedanken experiments. By this term, broadly speaking, one refers to rather ide- 
alized physical scenarios accompanied by reasonings aimed at highlighting some 
specific conceptual aspects abstracted from a difficult scientific problem, with 
disregard of any technological limitation. One of the most famous gedanken ex- 
periments was proposed by Schrodingcr to draw the attention to the paradoxical 
results to which a blind application of quantum mechanics to the macroworld 
may lead. This example is known as the Schrodingcr cat paradox and today it 
represents the archetype of all the physical situations wherein quantum super- 
positions of two macroscopically distinguishable states are taken into consider- 
ation. It is interesting to point out that some gedanken experiments, differently 
from the one involving the unlucky cat, are not only conceptual constructions 
but are also conceived as real, albeit impracticable, experimental schemes. 

The peculiar feature of a thought experiment of this kind is the striking 
contrast between the simplicity of the ingredients involved (a single atom, a 
single harmonic oscillator, a single quantized electromagnetic mode, etc.) and 
the unavailability of the necessary experimental setup. 

Fortunately, owing to technological progress, experiments considered unre- 
alistic up to recently, can now be carried out. As a result, various intriguing 
genuinely quantum effects traceable back to the superposition principle, entan- 
glement, quantum interference, etc., are now within experimental reach. Ex- 
periments of increasing difficultly in cavity quantum electrodynamics (CQED) 
over the last years, have, for example, made it possible to test fundamental 
radiation-matter interaction models involving single atoms. Moreover, cavity 
field states possessing remarkable nonclassical features have been generated and 
detected. 

In the domain of trapped ions, on the other hand, a single material quan- 
tum harmonic oscillator has been experimentally realized and the possibility of 
tailoring at will its coupling with electronic degrees of freedom has been practi- 
cally demonstrated. It is therefore to be expected that these two research areas 
will play a decisive role in resolving the basic paradoxical aspects of quantum 
physics. 

The scope of this paper is to review the state of the art of the theoretical 
investigations and of the experimental observations concerning the dynamical 
features of the coupling between single few-level atoms and two bosonic modes. 
The physical situations which we shall refer to, belong to the experimental 
domains of cavity quantum electrodynamics and of trapped ion physics. In 
the latter case we shall describe the physical properties of an ion harmonically 
confined in a bidimensional Paul microtrap. The motion of the center of mass 
of the ion is coupled to its internal dynamics using appropriate laser beam 
configurations. The electronic degrees of freedom of the ion are represented 
by pseudo-spin operators whereas the motional ones by bosonic creation and 
annihilation operators associated to the unperturbed bidimensional harmonic 
motion of the ion. 

The physical scenario relative to the problems we shall be faced with in a 
CQED context, involves one or more two- or three-level atoms interacting, one 
at a time, with two quantized electromagnetic modes sustained by a high-Q res- 
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onator. In this case too, the dynamical properties of the system arc investigated 
using hamiltonian models characterized by the presence of bosonic variables de- 
scribing the two quantized electromagnetic modes and of pseudo-spin atomic 
operators. Thus the two different research areas share, at least, the mathe- 
matical language and possibly, in addition, the potentiality for realizing, with 
the help of different experimental setups, classes of simple but basic pseudo 
spin-boson hamiltonian models. 

It turns out that this is indeed true, even if, as we shall see, the trapped 
ion context exhibits a greater practical feasibility which suggests experimental 
investigations beyond the possibilities of a typical CQED apparatus. 

What are the motivations which justify a review of these physical situations? 
In particular: what is the specific role of the two bosonic modes? 

A first answer is that the dynamics of several hamiltonian models describing 
such systems is exactly treatable and, in most cases, testable in the laboratory. 
A second more intriguing reason is that investigating these systems is likely to 
shed light on basic questions of quantum mechanics. The point to be appreci- 
ated is indeed that, studying such systems, one has the opportunity to induce 
entanglement and to control its evolution in a tripartite physical system atom - 
mode 1 -mode 2 with acceptable experimental efficiency. 

As we shall see, this entangled condition is essential for many theoretical 
proposals and experimental studies, both in CQED and trapped ion context, 
concerning the generation of a large class of interesting entangled states of the 
two bosonic modes only. We emphasize, in addition, that the presence of two 
modes in the physical systems under scrutiny, can be successfully exploited also 
to produce quantum target states of one mode only. In this last application the 
other mode plays an auxiliary role, leading to an effective improvement of the 
preparation protocol. 

The quantum states of the bimodal subsystem is intrinsically interesting be- 
cause, for example, it allows the production of bimodal states having nonclassical 
signatures. Of particular relevance is also the fact that the preparation of some 
entangled bimodal states yields the possibility of following the degradation of 
quantum coherence, under almost ideal conditions. 

The importance of experimental observations of this kind might be crucial to 
elucidate the nature of the decoherence mechanisms, unavoidably plaguing the 
time evolution of any mesoscopic or macroscopic physical system. Quite recently 
it has been possible to follow in the laboratory the degradation of mesoscopically 
sized Schrodinger cat like states advancing towards the exciting goal of realizing 
a gedanken experiment. 

We wish to emphasize that also the time evolution of the tripartite system 
is of interest in its own, at least under the aspect of bringing to light the mech- 
anism governing the onset and the time development of quantum mode-mode 
correlation effects. In this paper, for example, we shall review the so called 
parity effect recently discovered in the dynamics of a two-level atom quadrati- 
cally coupled to two bosonic modes. This intermode correlation effect is related 
to the granularity of the energy of a quantum harmonic oscillator. The oc- 
currence of some interesting consequences of the parity effect both in CQED 
and trapped ion contexts, shall be presented and discussed, clarifying in this 
way, the motivation for reviewing together phenomena belonging to these two 
different domains. 



This paper is organized in two large sections. The first is devoted to system- 
atically presenting the dynamical properties of the coupling between few-level 
atoms and two quantized electromagnetic modes. The second explores, in the 
same spirit, the physics of an ion confined in a bidimensional Paul trap and 
subjected to appropriate laser beam configurations. 



2 Atoms in bimodal cavities 



2.1 Introduction 

Over the last thirty years many theoretical investigations have been addressed 
toward the understanding of nonlinear effects occurring in the dynamics of an 
atom in a high-Q cavity. The interest toward this research area was mainly 
spurred by the large amount of experiments revealing the appearance of sub- 
poissonian photon statistics pi p[ a mong other intringuing features of quantum 
radiation-matter interaction |^f|4| g. Both theoretical and experimental ac- 
tivities have concentrated on trying to understand simple non trivial models 
of quantum optics involving a single atom, regarded as a few effective energy 
levels, and one or more near resonant modes of the quantized electromagnetic 
field. The prototype of such systems, proposed by Jaynes and Cummings in 
1963 |(|, describes a two-level atom resonantly interacting with a single- mode 
quantized field. When the Rotating Wave Approximation (RWA) is made, the 
model becomes exactly solvable and its dynamical features can be analytically 
brought to light revealing remarkable properties. A review of the Jaynes and 
Cummings model (JCM) has been presented by Shore and Knight in 1993 @. 
The discovery of interesting aspects of the JCM as well as the developments in 
CQED experiments involving single Rydberg atoms within single-mode cavities 
[^) , have stimulated an intense research devoted at highlighting and generalizing 
the original idea and physical scenario presented by Jaynes and Cummings. In 
particular their work has been extended to include more elaborate interactions, 
such as multiphoton || [To[ [il~[ o r intensity-dependent coupling |l2|, [l3|, |lj, [l5| , 
damping and dissipation |16|, |17| , multiple atoms |ll| [lj], , multilevel atoms 
[ [?l[ p2| and multimodc field |2j|, |4j. It is interesting to underline that histor- 
ically many extensions of the JCM have been reported as examples of exactly 
treatable nonlinear matter-field interaction models and that some of them, in- 
cluding the photon loss, provide realistic description of the experiments. 

After the experimental work of Brune et al ^5|, which have successfully 
realized a two-photon single mode maser operated in a high-Q cavity, a large 
amount of experimental and theoretical research has been done on the atomic 
two-photon transition processes in a cavity. The experiment successfully per- 
formed at the ENS laboratory, represents an important milestone in quantum 
optics because it has practically realized a two-photon cascade micromaser. Un- 
der the stimulus of the success of the experiment of Brune et al., other schemes 
for a two- mode two-photon maser have been presented P^| . 

More recently, much attention has been devoted to the generation of entan- 
gled states of two modes of two spatially separated cavities p7|| . Apart from 
their intrinsic theoretical interest, these states play an important role in the 
framework of the new field of quantum communication. For example, a method 
for generating a maximally entangled state of two cavities has been reported as 
an intermediate step in the teleportation procedure proposed by Davidovich et 
al [^8| . Recently, Zubairy et al |2{| have shown the possibility of teleporting a 
radiation field state from a cavity to another one, provided that the generation 
of a two-resonator entangled state for fixed number of photons inside the two 
cavities, is feasible. Notwithstanding the relevance of this topic, the study of 
protocols aimed at entangling two separate single mode cavities, goes beyond 
the scope of this topical review. 
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Section II of this paper is devoted to the interaction between bimodal cavities 
and few level atoms. In particular in Section II B a brief review of the JCM 
is presented whereas in Section II C its generalization to a bimodal interaction 
is discussed putting into evidence the interest stemming from the consideration 
of two modes instead of one. In Section II D a short discussion on two-photon 
micromaser is presented underlining its difference with respect to the one-photon 
micromaser. In Section II E interesting nonclassical dynamical effects of two- 
mode JC models are reported both in correspondence to the nondegenerate case 
(u>i W2) and in the degenerate one (u>i = UJ2) when one-photon processes are 
assumed. In Section II F we concentrate our attention on nonlinear interaction 
models for which the coupling mechanism between the atom and the bimodal 
cavity is based on two photon processes. In particular Section II F 1 is devoted 
to the so called Raman coupled model characterized by the fact that the net 
gain or loss of photons is zero in any process. The models adopted in Sections 
II F 2 -II F 4 are such that the respective basic processes involve the net gain 
or loss of a pair of photons having or not the same energy. Finally, in the last 
section (section II G) of this first part, several schemes aimed at generating 
nonclassical states in a bimodal high-Q microcavity are reported. 

2.2 A brief review of the JCM 

In the spirit of this topical review, we start devoting this section to a brief 
discussion of the JCM being it the simplest model describing the interaction 
between a single few-level atom and the quantized cavity field. This model has 
been extensively studied by many authors and numerous non-classical effects 
have been predicted and in some cases verified in laboratory [jlj 0, [|, ||, ||, ^, || . 
The JCM consists in a single two- level atom, whose two states are denoted by |+) 
and |— ), coupled to a quantized single mode of a cavity field having a frequency 
oj close to the transition frequency loq of the atom. The effective hamiltonian 
describing such a physical system, in the Rotating Wave Approximation, can 
be written in the form (h — 1): 

Hjc = i^oS z + loo) a + [Aa^S*- + h.c] (1) 

where a and denote the annihilation and creation operator of the cavity mode 
whereas the atomic degrees of fredom are described by the pseudospin operators 
S z and S± satisfying the following relations 

2&|±)=±|±>, S±|±) = 0, 5±| T ) = |±) (2) 

The coupling constant A depends on cavity volume, dipole strenghts and mode 
frequency |?], p0| . In the JCM the familiar Rabi oscillations characterizing the 
time evolution of the atomic inversion S z , are affected by the initial distribution 
of photon numbers |3l| . It has indeed been demonstrated that the interfer- 
ence between all the possible Rabi frequencies involved in the dynamics of the 
system, causes a dephasing or collapse of the Rabi oscillations. Moreover, the 
discrete nature of the photon number distribution leads to a partially rephas- 
ing or revival of such oscillations. The existence of this revival property has 
been experimentally observed and, being related to the discretness of the field 
excitations, provides direct evidence for the truly quantum nature of radiation 

@- 



2.3 Two-mode Jaynes Cummings models 



A natural generalization of the JCM consists in considering a single atom inter- 
acting with two modes of the quantized cavity field. When such a physical situa- 
tion is considered, several effective models can be constructed depending on the 
interaction mechanism characterizing the specific coupling. In particular, the 
Hamiltonian model describing the atom - two cavity modes coupling, assumes 
different forms depending on whether one-photon or two-photon processes are 
considered. In this last case the energy exchange between the atom and the 
bimodal cavity is mediated by virtual levels |32j not explicitely involved in the 
effective hamiltonian. In order to better clarify this point, consider one of the 
method followed by many authors for deriving effective models, namely the so- 
called "adiabatic elimination of virtual levels procedure" Q [34|, |3^, |3(], |37| . Its 
starting point consists in writing the full microscopic Hamiltonian in the dipolc 
approximation and neglecting all the counter-rotating contributions in accor- 
dance with the RWA. Then the Heisenberg equations of the field and atomic 
operators in the interaction picture are considered. If the population in a given 
atomic level is relatively much smaller than the other ones, then, formally inte- 
grating the relative Heisenberg equations, this level, called virtual atomic level 
[|D , can be adiabatically eliminated. 

It is important to underline that the elimination of these virtual or interme- 
diate states in the construction of two-photon effective models, leads to Stark 
shifts of atomic energy levels. As we will see, these shifts are proportional to 
photon numbers and act as intensity-dependent detuning. 

The interest toward two-mode JCM originates from the several remarkable 
features characterizing their dynamics. In our opinion one of the most interest- 
ing property of the two-mode JC models is that they provide the possibility to 
use one mode to modulate, amplify or control the other one. As we will show 
later, this property turns out to be very useful in the production of nonclassical 
states of a single cavity mode. Another relevant aspect characterizing such kind 
of systems lies in their wide potential tunability since oscillation or amplifica- 
tion can occur at any two frequencies u>\ and w%, with eventually w\ = ui2, such 
that u)\ + UJ2 matches the atomic transition frequency. Before concluding this 
section we wish to point out that in order to describe two-mode JC models, 
we will denote by and a\ [i = 1,2) the annihilation and creation operators 
respectively of the i— th cavity mode. Moreover, when the atom is effectively 
describable as a two- level atom the pseudospin operators S z and S±, such that 
the condition ^ are satisfied, will be used. 

2.4 Two-photon micromaser 

Micromasers are quantum oscillators operating on a Rydberg atomic transitions, 
namely corresponding to high principal quantum numbers. In a micromaser an 
atomic beam is injected in a high-Q microwave cavity. The atomic injection 
rate is low enough to ensure that there is just a single atom inside the cavity 
at one time. The experimental realization of these devices provides a very sim- 
ple way to test fundamental predictions related to radiation-matter interaction. 
In particular, when the basic interaction process between each atom and the 
cavity field involves the emission or absorbtion of one photon per atomic tran- 
sition (one-atom micromaser) , the micromaser provides a near-ideal realization 
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of the problem firstly theoretically treated by Jaynes and Cummings. These de- 
vices also allow the operation of two-photon oscillators involving a degenerate 
or nondegenerate two photon transitions between two atomic levels of the same 
parity. These two-photon micromasers present additional interesting features 
with respect to the one-photon micromaser, which are connected to the starting 
up of the system. Classically, the zero-field state is always stable so that an 
injected field is required in order to build up the oscillation. However, quantum 
fluctuations may drive the occurrence of a vacuum state instability, thus in- 
ducing the system to begin the oscillation. The two-photon micromaser indeed 
spontaneously starts without any triggering required. Moreover the two-photon 
micromaser provides a clean setup for studying the two photon processes. A 
two-photon micromaser was realized by Brune et al. []38f using a niobium super- 
conducting cavity with Q ~ 10 s , resonant at u> ~ 68GHz. They in particular 
exploited a two-photon degenerate transition in Rubidium atoms. 

When the two photons are emitted into different modes of the cavity (nonde- 
generate two-photon transition) important features of the two-photon process, 
such as for example the correlation between the two photons emitted in two 
different modes, can be studied. 

It is important to emphasize that the mere fact that the atomic system emits 
or absorbs photon pairs leads to a field statistics very different from the one of 
the ordinary laser or maser fields and makes it a privileged candidate for the 
generation of squeezed light. In such systems the Rydberg atom plays a role 
similar to a coupler which correlates the two modes of light changing the photon 
statistics inside the cavity. 

The main novelty of the nondegenerate two-photon micromaser (two differ- 
ent cavity modes of frequencies u>\ and uji respectively with uii 7^ CJ2), with 
respect to the degenerate one (lu\ = w 2 ), is the strong correlation between the 
two different modes, which may produce a 50 % squeezing in the difference of 
intensities at steady state p6fl . 

2.5 One-photon processes 

In this section we consider the coupling between few atoms and two modes of 
a high-Q cavity field assuming that the interaction mechanism relies on single- 
photon processes. In particular in the subsection II E I we focus on the case 
of an effective three-level atom and two nondegenerate modes of a single cavity 
whereas in subsection II E 2 we assume that the two cavity modes are degenerate 
and that the atom can be effectively represented as a two-level system. 

2.5.1 Three-level atom interacting with a bimodal cavity 

A natural extension of the JCM consists of a three-level atom coupled to one 
or two modes of a high-Q cavity field. This generalization allows, for example, 
the investigation of physical phenomena associated with two-photon processes. 
A general treatment of the theory concerning such systems has been presented 
in |40| where arbitrary detunings are considered. Denoting the three atomic 
levels by |0), |I) and |2) and limiting our discussion to the electronic dipole 
transition |0)-|1) and |1)-|2), there are three possible distinct atomic level con- 
figurations as shown in fig. [|. Following fljlf , these configurations are called A, 
S and V . In this context the transition |0)-|2) is forbidden. 



When the interaction between the atomic system and the two cavity modes is 
considered, it is interesting to study the effects of intermode field correlations on 
the dynamics of the atomic system and viceversa. This particular aspect of the 
problem forms the basis of the paper of Lai et al |Q . In order to investigate on 
the occurrence of dynamical effects related to the presence of initial correlation 
in the physical system, the authors consider an SU(2) coherent state 43 as 
an example of a two-mode correlated state and a three-level atom in the A 
configuration as atomic transition sequence. The SU(2) coherent state has the 
property that, indicating by m the maximum possible number of excitations 
in each mode, if n photons are observed in one mode then (m — n) photons 
will be observed in the other one. This means that the SU(2) coherent state 
is an eigenstate of the total excitation number operator a\a\ + a\a,2 relative to 
the eigenvalue m. The three-level atomic system consists of two dipole allowed 
transitions, |0)-|1) and |1)-|2) respectively (see fig. 0(a)), resonant with two 
different modes of the quantized radiation field. The interaction Hamiltonian 
describing the coupled system, in the RWA, can be written down as 

= 5 i(a t 1 |0)(l| + h.c.) + ff2 (4|2)(l| + h.c.) (3) 

where g\ and g% are the coupling constants relative to the two different channels 
|0) — |1) and 1 2) — |1) respectively. The particular form of initial field correlation 
imposed by the authors to the two cavity modes, leads to generalized two- 
photon Rabi frequencies that are independent of the total photon number m, 
provided that the atom-field coupling constants are set equal. Thus, under this 
condition, the atomic population probability exhibits purely periodic behaviour. 
Moreover the revival times, that is the time instants at which the oscillations 
revive, are much shorter than the ones exhibited assuming the field in the two- 
mode uncorrelated binomial state. This behaviour is of relevance in view of the 
fact that the photon statistics of the individual modes are the same in the two 
cases considered by the authors. 

The A configuration is also exploited in ref. p3j . In this paper the authors 
define the slowly varying complex amplitudes of the two- mode field as following: 

rf'i4E (% e<W3 " + aK*"*) (% e<W3 " - a}e-^*) (4) 

22 3=1,2 22 3 = 1,2 

where u)j is the frequency of the j— th mode. By definition [e^cfe] = \ so that 
squeezing of the two-mode field occurs whenever one of the amplitude operators 
satisfies the condition 

< (Arf,) 2 >< \ (5) 

In particular the authors demonstrate the occurrence of two-mode squeezing 
starting from the two cavity modes prepared in a coherent state. Other inter- 
esting phenomena are found as, for example, the periodic recovery of squeezed 
states caused by the involved two-photon processes. 

The interaction between a three-level atom and two different cavity modes 
has been shown to produce also normal as well as higher-order squeezing. In ref. 
[fl5|, M.H. Mahran studies the interaction of a three-level V-shaped atom (see 
fig. |l|(c)) with two quantized electromagnetic modes discussing the possibility 
of occurrence of squeezing. Following the definition proposed by Loudon and 
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Knight Eq] , the author has calculated the two-mode squeezing when the atom is 
initially prepared in its ground state. The results obtained by Mahran brought 
to light the dependence of the two-mode squeezing on the initial mean photon 
number. In particular the observed amount of squeezing decreases as the mean 
value of photons initially present in the two modes increases. Finally the author 
considers the sum squeezing, as defined by Hillery j47j, for an atom initially 
prepared in its ground state, proving that such a quantity is smaller than the 
amplitude-squared squeezing. 

The ladder or S configuration, displayed in fig. 0(b) is investigated in jl8| 
bringing to light that it can produce very good quantum nondemolition mea- 
surement performances. 

Phase properties of the field in the two-mode three-level problem are in- 
vestigated inj49|]. In particular, using the Hermitian phase formalism of Pegg 
and Barnett [pOj, phase properties of each individual mode as well as their joint 
probability distribution and correlation function are studied. The authors show 
that phase properties of the field reflect the collapses and revivals of the atomic 
level occupation probabilities except for the Raman scattering in correspondence 
of which the collapses and revival of the atomic population can be completely 
decorrelated from the phase of the field. 

Very recently Q the master equation for the reduced density operator of the 
field in the two-mode micromaser pumped by ultracold A-type three-level atoms 
has been derived and the steady state photon probability distribution under the 
principle of detailed balance has been evaluated. Analogously to the definition 
given by Scully et al. the system now described is called two-mode mazer 
(microwave amplification by the z motion induced emission of radiation). When 
indeed the micromaser is pumped by laser cooled atoms, quantization of atomic 
external motion becomes necessary. The interesting result obtained by the au- 
thors is that the degree of anticorrelation between the cavity modes increases 
when the micromaser is pumped by ultracold atoms instead of fast atoms. This 
leads to much stronger sub-Poissonian statistics for each mode. 



2.5.2 Degenerate one-photon two-mode model 

In this section we focus our attention on the one-photon interaction between a 
single two-level atom and two degenerate modes of a lossless cavity described 



with the help of the following Hamiltonian model |53 



H d = u V a% + u) S z + V g^e-^&l S_ + h.c.) (6) 



fj,= l,2 /-i=l, 2 



This model conserves the total excitation number N — a\a\ + a 2 a>2 + S z + ^. 

2 , 

Moreover it possesses the purely bosonic constant of motion C = g -i+ g -i a\a\ + 

2 gl i atci2 9l , 92 2 (olc^e'^ 2 "^ 1 ) + h.c). The existence of these two constants 

91+92 2 Bi+gi v 1 ; 

of motion leads to the possibility of unitarily transforming into a system 
consisting of an one mode JC subsystem and a decoupled other new bosonic 
mode. Thus the dynamics of the model given by equation (^|), may be exactly 
solved starting from an arbitrary initial condition. 

When one mode, say the first, is prepared in a Fock state with n ^> 1 
excitations and the other one is empty, the field dynamics is dominated by a set 
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of n incommensurate Rabi frequencies leading to the occurrence of oscillatory 
net exchanges of a large number of photons between the two modes. These 
oscillations of the field populations display amplitude decay and, over a longer 
time scale, exhibit revivals and collapses. This behaviour is displayed in fig. |^ 
where we have in particular chosen n — 20. Analogous phenomena take place 
in the dynamics of the system when the initially populated mode is prepared 
in a coherent state [Q. The reaction of the atomic dipole on the field and the 
granularity of the field itself, are at the physical origin of the highly nonclassical 
regime attained by the system under scrutiny. In J54j the role played by the 
atom-cavity field entanglement is brought to light and related to the occurrence 
of a cooperative intermode effect which is shown to be incompatible with a 
neoclassical approach. 



2.6 Nonlinear interaction models 

Over the last two decades the interest toward the dynamics of nonlinear models 
effectively describing the interaction between an atom and two quantized modes 
of a high-Q cavity, has grown up in a significative way. With the emergence 
of the two- mode two-photon micromaser of Laughlin and Swain |]55| , it has 
become possible to experimentally observe certain nonlinear effects including 
field squeezing and emission spectra. Both theoretical and experimental activi- 
ties concentrate on trying to understand simple non trivial models of quantum 
optics involving single effective two-level atoms and two independent modes 
of the radiation field, having or not the same frequency, coupled by means of 
two-photon processes. In this section we focus our attention on nonlinear two- 
mode models, bringing to light the main features characterizing their relative 
dynamics. 



2.6.1 Raman coupled model 

Let's consider two quantized cavity modes of frequency uii and u>2(< wi) respec- 
tively and a three-level Rydberg atom whose energy diagram is given by fig. [|. 
Following Gerry and Eberly the level 2 can be adiabatically eliminated 
provided that the condition of large detuning A ^> E 3 — Ei is satisfied. 

Under this hypothesis the effective hamiltonian describing the system under 
scrutiny can be written down in the following form: 

H R = + H? (7) 

with 

Hq = uj S z + ujia\ai + uj 2 'a\^ Hf = gRCiia\S + + h.c. (8) 

The resulting model therefore consists of two-nondegenerate "ground" states 
connected by zero-photon processes where a photon absorbed in one mode is 
emitted in the other one with the atom making transitions through the virtual 
state |2)[|5). 

The system under scrutiny, effectively described by Hr, may be interpreted 
as a cavity version of Raman scattering in which mode 1 is the pump field, mode 
2 is the Stokes field and the anti-Stokes mode is eliminated by setting a cavity 
off-resonance condition |36). As a direct consequence of the fact that we are 
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dealing with zero photon process, the total photon number operator defined as 
N = a\ai + a\a,2 is a constant of motion J35|. 

The two-mode Raman coupled model considered in this section, has been ex- 
tensively investigated showing the existence of collapses and revival phenomena 
[ p5| , resonance fluorescence |57| , population trapping K(| and several nonclassi- 
cal properties of the electromagnetic field such as the appearing of antibunching 
light, violation of the Cauchy-Schwartz inequality and squeezing |35|, ^8|. As far 
as the collapse and revival phenomena of the Rabi oscillations, it has been proved 
that, differently from the JC model, in this case it can be periodic |55| (58|. In 
particular, assuming coherent states for both modes, the atomic inversion be- 
tween the two nondegenerate ground states exhibits periodic revivals but with 
many secondary revivals as well. These last revivals become less significant for 
higher initial average photon numbers. The physical origin of such collapse and 
revival patterns was partially brought to light in Q but it has been more stud- 
ied, at least for the case of equal average photon numbers in both modes, by 
Cardimona et al. [[37] . As far as nonclassical properties of the electromagnetic 
field are concerned, Mahran [^9), has, for example, studied the production of 
squeezing. 

The Hamiltonian model (^)- (||) adopted in |35, 37, 56[ [59| neglects the Stark 



shift terms arising from the adiabatic elimination of a virtual level. Relaxing this 
approximation implies indeed that the resulting two-level system is no longer 
two-photon resonant but it is detuned by an intensity-dependent amount. In 
other words the unperturbed atomic hamiltonian u)qS z appearing in , has to 
be substituted by 

(E 3 + |4« 2 )|3)(3| + (E 1 + ^a\*i)m\ (9) 

where <?j (i = 1, 2) is the atom i — th mode coupling constant. Buzek and Knight 
[ p0[ have shown that this substitution may lead to quite different evolutions 
even for low-intensity fields. The effects of the Stark shift terms on the atomic 
inversion dynamics when a multilevel atom interacts with a two-mode quantized 
radiation field have been also investigated by Cardimona et al. j6lj. They have 
shown that these Stark terms are relatively unimportant for Raman scattering 
but could not be ignored in equal- or unequal-frequency absorption. 

Very recently a class of coupled-channel cavity QED models describing two- 
photon resonant Raman processes in multiple- A configuration has been solved 
analytically by an algebraic method based on the introduction of an angular 
momentum algebra for this multiwave mixing system ||62| . 

2.6.2 Nondegenerate two-mode two-photon JC model 

The system we consider in this subsection is a generalized JC model in which 
the atomic transitions are mediated by nondegenerate two-photon absorption or 
emission. As it is well known, in the two-photon processes at least an interme- 
diate level, assumed to be coupled to the ground and excited states respectively 
by dipole allowed transitions, is involved. As shown in Fig. |] we denote by E g , 
Ei and E e the lower, intermediate and upper atomic energy levels and, indi- 
cating by u>i and 0J2 the frequencies of the two cavity modes, suppose that the 
conditions Ei — E g = uj\ — A, E e — Ei = L02 + A are satisfied. It is possible 
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to demonstrate that, if |A| is much larger than the one-photon Rabi frequency 
of the oscillations between \i) and |<7) and between \i) and |e) respectively, the 
effective Hamiltonian model describing the system under scrutiny assumes the 
form: 

H nd = uia\ai + Lo 2 a\a 2 + {u>o + P 2 a\a 2 - (3\a\a\)S z + (Xaia 2 S + + h.c.) (10) 

The alkali-metal Rydberg energy diagram, for example, provides concrete 
examples fulfilling the conditions leading to an atom - bimodal cavity interac- 
tion Hamiltonian such as H n d |25|, Pqj . Observe that the hamiltonian model 
given by equation (|Io|), contains an intensity dependent detuning described 
by the two terms proportional to f3\ and (3 2 respectively, whereas the constant 
A measures the effective two-photon atom-field coupling. This Hamiltonian is 
exactly solvable and a unitary operator accomplishing the canonical dressing of 
the two-level atom by the bimodal cavity field, has also been constructed 

The Hamiltonian model H n d was firstly introduced by Gou |64|, |65| . In 
neglecting the Stark shift arising from the presence of virtual levels, that is 
putting in equation ( |Io| ) 0i ~ (32 ~ 0, the quantum dynamics of the atom 
and the temporal evolution of the field statistics were studied. In particular 
the author considers two types of initial state of the field, namely ordinary 
uncorrelated coherent states in both modes and the two-mode squeezed vacuum 
state where the number of photons in one of the modes is equal to the number 
of photons in the other one. It has been noted that the dynamics of the model 
is qualitatively different in correspondence of these two kinds of states. In more 
detail it has been shown that new features of quantum collapse and revival 
occur in presence of the two-mode squeezed vacuum. The revivals here appear 
periodically sharp. The origin of such a behaviour is strictly related to the 
particular class of Rabi frequencies involved in the dynamics of the respective 
physical systems. If, on the contrary, the field is initially prepared in a product 
of two coherent states, the atom exhibits the phenomena of collapse and revival 
as well, but no periodicity can be observed. The different behaviour of the 
atomic dynamics in correspondence to initial correlated and uncorrelated cavity 
field has also been put into evidence by Joshi and Puri |66| as well as by Gerry 
and Welch [p7| . In both papers the Stark shift terms are neglected. 

In |68| Ashraf has investigated the cavity field spectra for various initial fields 
of a nondegenerate two-photon JC model neglecting the two terms proportional 
to Pi and (3 2 respectively. 

The results presented in |54| [37[ |6|] could be qualitatively or quantitatively 
affected by the presence of the Stark shift of the two effective two atomic levels. 
As shown by Gou [[35), indeed, in most circumstances these terms cannot be 
ignored, even if for certain correlated two-mode field states, the elimination 
of Stark shifts is possible at least when one evaluates the ensemble averages 
relevant to the phase insensitive terms. However, for the phase sensitive averages 
such as the squeezing of cavity fields, the Stark shifts still cannot be ignored 
in any case. Also Cardimona et al |3l] have proved that, differently from the 
Raman coupled model, for unequal-frequency absorption or emission, the Stark 
effect could not be negligible. 
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2.6.3 Intensity dependent two-mode two-photon JC model 



In this section we investigate on the nondegenerate two-mode two-photon JC 
hamiltonian model assuming that both the coupling strcnght and the atom-field 
detuning may depend on the populations of the two modes. The hamiltonian 
describing such a general nonlinear model may be written down as J63[: 



where the RWA has been used. The two functions G{a\a\, a 2 a,2) and F(a\ai, a\a<z) 
appearing in equation (|l2|), are introduced to take into account the pos- 
sibility of intensity dependence both in the atom-field detuning and in the 
coupling constant. Let's observe that putting in particular u> G({alai}) = 
uiq — (3ia\ai+ (32^20,2 and F({a\ai}) — I we recover the nondegenerate two-mode 
two-photon JC model discussed in the previous section. In ||63f a symmetry- 
based detailed construction of the atom-cavity eigenstates as well as of the 
energy spectrum of the dressed atom is presented. The starting point of such 
a construction consists in finding canonical symmetry transformations for Hid 
and the related constants of motion. In particular they prove that the opera- 
tors N = a\a\ + a\a,2 + 2S Z + 1 and D = a\ai — a\a,2 are constants of motion 
indipendently from the specific form of the two intensity dependent functions 
F and G. The existence of these two constants of motion turns out to be very 
useful in the construction of a unitary operator diagonalizing Hid- It is worth 
noting that the symmetry-based conclusions envisaged in p| may be applied 
whatever the analytic dependence of the intensity-dependent terms on the pop- 
ulations of the two-cavity modes is. This result allows to highlight the role of 
the nonlinear contributions on the two-level atom dressed by the cavity field. In 
order to do this, in p3] the authors study the quantum dynamics of the initially 
excited two-level atom interacting, in the context of the general model, with a 
cavity field prepared in the tensorial product of a coherent state and a Fock 
state. Exploiting the flexibility of their results the authors compare the atomic 
evolution of several particular models finding a rich dynamical behaviour dom- 
inated by the presence of collapses and revivals of the Rabi oscillations. The 
occurrence of these collapses and revivals is a common phenomenon to all the 
models considered in | ]63| even if some other distinctive features characterize the 
time evolution of the atomic inversion in each individual model. 

The existence of sensitivity to the population of the field mode initially 
prepared in a Fock state is also pointed out. This effect is of quantum nature 
and it is particularly noteworthy since, as showed in the subsequent subsection 
II G 2 it leads to proposals for generating number states in a single cavity mode 
starting from the vacuum. 

2.6.4 Degenerate two-mode two-photon JC model. Parity Effect 

Suppose now that the two modes of the cavity field have the same frequency 
lo assuming that they differ either by the polarization vectors or by the direc- 
tion of propagation. The degeneration condition imposed to the two quantized 



Hid — H i + H 2 + i?3 



(11) 



with 



Hi = LUid\ai+LU2ala2] Fl 2 
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electromagnetic modes makes the effective hamiltonian model describing the 
two-photon interaction between the bimodal cavity and a single few level atom, 
very different from that adopted in the previous section. The effective model 
obtained in this degenerate case is not a trivial generalization of the nondegen- 
erate one ([lO]). The explicit construction of the effective Hamiltonian for the 
degenerate quadratic atom-held interaction leads to the following model: 

H d = H d + H d R + H d s + Hf (13) 

with 

Hq = uj S z +u! ^L^; Hr = (ria 1 al+r 2 a 1 alS z )+h.c; Hg = sS z ^ a 

(14) 

and 

&i = ( S X & + 9^2)S+ + h.c. (15) 

fi=l,2 

This Hamiltonian model has been recently introduced by A. Napoli and A. 
Messina |69| . The structure of Hd reflects the fact that the atom- mode- mode 
field energy exchanges in the degenerate case may take place following all the dif- 
ferent possible channels compatible with the conservation of the total excitation 
number operator = d\ai + a 2 a,2 + 2S Z + 1. It is easy,indeed, to demonstrate 
that the operator N commutes with Hd whatever the coupling constant s, r\, 
T"2, Ai, A2 and g are. The term of Hd depending on the constant g, for example, 
arises from the possibility that the atom makes transitions involving a simulta- 
neous gain or loss of one photon in both modes. The contributions related to the 
coupling constants T\ and r-i describe, on the other hand, mode-mode energy 
exchanges not involving the atom or changes in its state. These terms are known 
as Rayleigh terms. Stark terms describing a field-dependent detuning appear 
in Hd with the coupling constant s. Finally the possibility of two-photon ex- 
changes between the atom and each cavity mode is described by the terms of Hd 
containing a„/S+. The strength of this specific effective atom-field interaction 
mechanism is measured by A^. 

In what follows we concentrate our attention on the case correspondent to 
g = and Ai = — A2 = A, that is on all the physical systems effectively describ- 
able by means of the following Hamiltonian model: 

H d = uj S z +uj ^2 a^a^ + sS z ^ a^a^ + \(na-ia\ + r 2 aia\S z ) + ft(dp) 

H=l,2 A*=l,2 

+ [A(a x — a 2 ! )aS' + + h.c] 

In fjg| it has been shown that specific physical situations in correspondence of 
which the conditions g = and Ai = — A2 = A are satisfied, may be envisaged. 

It is possible to convince oneself that in this case the Hamiltonian Hd pos- 
sesses another hermitian constant of motion, namely C — a\a,2 + aia\. The 
existence of the two constants of motion, N and C, is at the basis of the pos- 
sibility of exactly evaluating the dynamics of the system bringing to light the 
existence of an intrinsically quantum phenomenon, christened in |69f| "parity 
effect" , strictly connected with the granularity of the radiation field. In order to 
describe this nonclassical effect, let's suppose to inject a two-level atom, initially 



prepared in its ground state, into the cavity where one mode is excited in a state 
exactly containing n photons whereas the other one is left in its vacuum state. 
After an appropriate n— dependent interval of time, the system described by Hd 
exhibits a peculiar sensitivity to the parity of the initial total number of pho- 
tons n. Stated another way, the physical system distinguishes between the two 
different initial conditions n odd or n even, showing macroscopically different 
quantum behaviours. This sensitivity to the parity of the initial photon number 
clearly manifests itself in the dynamics of the mean value of the population op- 
erators a^a^. In Fig. || we in particular display the time evolution of the /i = 1 
mode average photon number, correspondent to the initial conditions n = 20 
(Fig. ||(a)) and n = 21 (Fig. ||(b)) respectively, putting r\ = r 2 = s = 0. These 
plots show that the initial behaviour of the temporal variation of the photon 
number in the mode \i = 1 is characterized by a rapid reduction to a half and 
it is followed by a relatively longer interval of time dominated by a negligible 
photon exchange between the two modes. After such a "lethargy period" , the 
time evolution of the mode population appears to be very sensitive to the parity 
of the initial photon number. Indeed, choosing the initial condition n = 20 the 
field mode completes the transfer of its population to the other mode provoking 
an inversion in the dynamical role between the two modes. On the other hand, 
the plot of the mean value fhi{i) of the operator fii — a\ai, relative to the ini- 
tial condition n = 21, shows that, after the "lethargy period", the field mode 
reabsorbs the photons already transferred to the initially empty mode. 

The existence of this interesting nonclassical effect, that is the sensitivity to 
the parity of the initial number of photons present in the excited field mode, 
is a peculiar feature of the atom-degenerate bimodal cavity field interaction 
described by Hd and, as demonstrated in J69[ , it does not appear related to the 
simplifying assumption r% = r-i = s = 0. 

We wish to underline that, as we will show in the next section, the parity 
effect can be successfully exploited for generating highly nonclassical bimodal 



Before concluding we wish to anticipate that in the third section of this 
paper, the parity effect shall appear also in the dynamics of some hamiltonian 
models of interest in trapped ion context. Such an occurrence is by no means 
surprising in view of the fact that, at a formal level, CQED and trapped ion 
domains make use of similar models. The important point to stress is however 
that, due to the different meaning attributed to the same hamiltonian model 
within the two experimental realms, very different dynamical consequences of 
the parity effect in these two context, shall be brought to light. 

2.7 Generation of nonclassical field states in a bimodal 
high-Q microcavity 

Investigating the dynamical behaviour of a physical system from a theoretical 
point of view, very often involves more or less tacit assumptions about the pos- 
sibility of preparing the system in a given state at t = 0. In this sense the 
experimental realization of assigned initial conditions for a system is a prereq- 
uisite for studying fundamental aspects of quantum theory ||, Q . Besides this 
inherent theoretical interest, the state preparation plays a central role also in the 
context of the new field of Quantum Information theory encompassing Quantum 
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Computing and Quantum Cryptography |73, [74], [75], f76[ . 

So far two approaches have been mainly exploited in order to produce desired 
prefixed states. The first one consists in finding the appropriate Hamiltonian 
which can transform a given initial state into a desired one via unitary evolution. 
The other commonly used approach consists in making a quantum measurement 
on one part of an entangled system provoking the wave function collapse of the 
other subsystem. 

Over the past few years, also under the stimulus of this rapid development in 
the experimental techniques, several ingenious procedures for generating non- 
classical states of the electromagnetic field have been reported in the realm of 
cavity quantum electrodynamics [Q, |77|, fr| @ fc], || |f , f| || ||, f(| §7|, || . 
Particular attention has been devoted to the creation of single-mode number 
states in view of the fact that they are considered as prototype of nonclassical 
state j77], |7?| |8(], ||] ||[ ^5], |86|. These states, corresponding to an exact num- 
ber of photons, possess in addition maximum information capacity and therefore 
they might be used for encoding and processing quantum information. 

A lot of attention has moreover been devoted to the possibility of preparing 
superpositions of two or more macroscopically distinguishable states j70], Jn], [78], 
Q |8£|. Such states, called Schrodinger cat like states, provide an useful starting 
point to progress with the understanding of decoherence phenomena and to gain 
more insight about the elusive border between classical and quantum worlds Q . 

Most of the proposals present in literature deal with the generation of non- 
classical states of a single cavity mode [Q, |77[ [7^, |o[ ||l], |83|. The interesting 
dynamical properties manifested by the two-mode models discussed in the pre- 
vious sections, suggest that also the generation of two-mode field states might 
play an important role in highlighting fundamental aspects of quantum theory 
In the present section we will show that the interaction between single atoms 
and two independent modes of a bimodal cavity can be successfully exploited 
to generate several nonclassical states of the bimodal field. At the same time, 
resuming an idea already present in the paper of Gou , we will demonstrate 
that the bimodal cavities can be also used for producing, with reasonable simple 
protocols, nonclassical states of a single mode of the cavity. In particular, in 
subsection II G 2 we will prove that when the attention is focused on a single ra- 
diation mode, the presence of the other one plays the relevant role of a "control 
subsystem" . 

2.7.1 Preparation of bimodal field states 

The basic idea of most of the experimental procedures aimed at generating 
assigned nonclassical states of one or more modes of a single cavity field, can 
be traced back to the possibility of manipulating the statistical properties of 
the radiation field by injecting, one at time, a flux of Rydberg atoms into the 
cavity. 

In the framework of micromaser-based preparation of cavity field states, 
one of the methods generally used for guiding the field toward a prefixed state 
is the so-called conditional measurement approach which consists in selecting 
only those sequences wherein each atom exits the cavity in a given state and 
discarding all others 78, [79], |o). Although most of the experimental schemes 



exploiting this idea are conceptually very attractive, nevertheless they are often 
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not mature enough to be qualified as realistic projects. Generally speaking, 
indeed, the high degree of precision assumed for the experimental apparatus 
as well as the fact that, even in such idealized conditions, the probability of 
success of the experiment is often very low, seriously counter the possibility of 
practically implementing these proposals. 



a) Correlated bimodal states from atomic coherences 

Consider, for example, the controllable generation procedure aimed at pro- 
ducing a class of arbitrary two-mode field states in a cavity, recently proposed by 
Deb et al |H| . Their proposal is based on the possibility of transferring atomic 
coherences to the cavity field exploiting the Raman coupled model. The au- 
thors show that, injecting into the cavity N three-level atoms each one initially 
prepared in a suitable superposition of its two nondegenerate ground states and 
measuring their internal state after the interaction with the cavity field, it is in 
principle possible to guide the field toward a state \<f>) = Yl n =0 7T n\ n > N ~ n ) 
starting from \ip(0)) = \0,N). Of course \<p) is normalized and its probability 
amplitudes is determined by the coupling mechanism and by details of the 
protocol. 

The idea firstly exploited by Vogel et al. |9^] of transferring atomic co- 
herences to the field, has also been used by Zheng and Guo [^3| who propose 
a rather formal way of preparing multimode field states having equal photon 
numbers in these modes. 

The nonresonant interaction of a two-level atom and a bimodal cavity field is 
exploited by Guo and Zheng [Q to propose a method for generating, at least in 
principle, any linear superposition of product of coherent states of the bimodal 
field. The basic idea of their scheme is that the photon statistical distribution 
of the field will not be changed by the interaction if the detuning between the 
atomic transition frequency and the frequency of each mode is much larger than 
the coupling constant. On the contrary the phase of the field evolves with the 
time. The authors claim that, by state-selective measurements on the atom, it 
is possible to produce entangled coherent states of the bimodal field provided 
that both the two modes are initially prepared in a coherent state. 



b) Pair Fock states 

In P4I the theory underlying a new experimental conditional measurement 
scheme for the preparation of selected equal-intensity bimodal Fock states is 
presented, carefully discussing its practical feasibility. The existence of a reli- 
able experimental method to produce pair-Fock states of controllable intensity 
besides its inherent theoretical interest, is also of relevance from an applicative 
point of view. In the field of quantum communications, in fact, it has been very 
recently analyzed the possibility of transmitting, with relatively high capacity, 
information by multiple parallel "number state single channels" [^5| . The ex- 
perimental project presented in [ p2| is based on the nondegenerate two-mode 
two-photon JC coupling, described by equation (|l0|), between single Rydberg 
atoms and a bimodal cavity. A peculiar and noteworthy characteristic of this 
specific atom-two mode field interaction mechanism is the possibility of obtain- 
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ing, starting from an empty cavity, pair-Fock states of the bimodal field with a 
relatively high intensity. The building up of an equal-intensity pair-Fock state 
of the microresonator field, in accordance with the proposal under scrutiny, re- 
sults from the success of a multiple iterative event. Suppose the cavity prepared 
in its vacuum state |0, 0) . A single excited atom enters the cavity where it 
simultaneously exchanges one photon with each mode as effectively described 
by equation (|l0|). Immediately after the atom has left the cavity its internal 
state is measured using an appropriate detection technique. It is possible to 
convince oneself that an atomic ground state outcome disentangles the atom 
from the field which, in turn, collapses in the state |1, 1). At this point a second 
atom is injected into the cavity with an appropriate velocity. If even the second 
atom is detected in its ground state exiting the cavity, then the field in the mi- 
croresonator is projected onto its Fock state |2, 2). It is not difficult to convince 
oneself that this procedure may be iterated so that a successful sequence of N 
atoms detected in their ground states immediately after they have left the cavity, 
builds up the equal- intensity pair-Fock state \N,N) of the bimodal field. The 
authors present a detailed quantitative analysis of the scheme showing, indeed, 
that its probability of success is of practical interest also in correspondence to 
high intensity pair Fock states (N ~ 20). On the other hand, notwithstanding 
the lowering of such a probability due to the presence of important technological 
limits of the experimental apparatus, the scheme maintains a reasonable com- 
patibility with the presence of sources of imprecision in the concrete realization 
of the experiment. 

A common aspect to many micromaser-based proposals for generating non- 
classical states of the cavity field is that they require the availability of ideal 
devices with which the interaction time between each atom and the field may 
be sharply selected. In connection with the generation of bimodal Fock states, 
the presence of unavoidable fluctuations in the atom-cavity interaction times is 
taken into account from the beginning in j85|. The theory developed in this 
paper clearly shows that the impossibility of sharply controlling the interaction 
time between each atom and the bimodal cavity, as well as the presence of other 
important sources of experimental imprecision, reduce the probability of success 
for the preparation of equal-intensity pair Fock states to values that, however, 
are still of experimental interest. As stressed in |8{|, this circumstance is strictly 
related to the interaction mechanism chosen for manipulating the cavity field. 



c) Maximally entangled states 

The nondegenerate two-mode two-photon JC coupling given by equation 
( |l0| ) appears to be a good basis for generating other classes of purely quantum 
states as, for example, maximally entangled two-mode states having the form 
\cf>±) = -^t|[|0, 0)±e 4 ^|l, 1)] |)(|. An attractive and relevant aspect of the scheme 
reported in pq] , is that the phase factor <fi can be manipulated from — 7r to 7r 
by simply varying an experimentally easily controllable parameter. 

Very recently, the possibility of entangling two independent cavity modes 
has been experimentally demonstrated [p8| . In particular, in |8£| the authors 
show that, using a single circular Rydberg atom, two quantized modes of a 
superconducting cavity can be prepared in the following maximally entangled 
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state characterized by the fact that the two modes share a single photon: 



-L( e ^0,l> + |l,0» (17) 

This EPR pair state is then revealed by a second atom probing, after a delay, 
the correlations between the two cavity modes. 



d) Photon number statistics invariant manipulation 

Exploiting the same interaction mechanism envisaged in the previous section, 
it is also possible to modify only the phases of the initial probability amplitudes 
of finding a well-defined population in the bimodal field simply injecting a single 
two-level Rydberg atom into the cavity p\ . As shown in pTj, this kind of 
radiation field manipulation can be useful, for example, in the construction of 
states orthogonal or quasi-orthogonal to the initial one or for controlling the 
squeezing parameter of a given single-mode state. 



e ) Schrddinger cat states from parity effect 

The existence of the intrinsically nonclassical "parity effect" discussed in sec- 
tion II F 4, is at basis of the possibility of generating entangled bimodal states 
having the nature of Schrodinger cats 70 , 7l| . We wish to underline that in this 



case the two modes of the bimodal cavity are degenerate and the hamiltonian 
model describing the two-photon interaction between the field and a single two- 
level Rydberg atom is given by equation (|l6|). Suppose to prepare the cavity 
leaving one mode in its vacuum state and exciting the other one in a coherent 
state | a). Such an initial condition can be easily realized in laboratory. The 
vacuum state is, indeed, reached controlling the temperature at which the ex- 
periment is performed 0. On the other hand, it is possible to excite one mode 
in a coherent state having a prefixed intensity, using classical currents Q. 
The experimental protocol described in fn], starts injecting in the resonator the 
two-level atom prepared in its ground state. In order to elucidate in which way 
the parity effect previously recalled may be exploited for reaching the prefixed 
scope, it is convenient to cast the initial condition imposed to the system as 
sum of the following two contributions: 

|^(0)) = e~^- Y -p=|2n, 0, -) + e"^" y |2n + 1, 0, -) 

(18) 

In view of the linear character of quantum mechanics and as a consequence 
of the parity effect, the initial state evolves toward two macroscopically differ- 
ent states respectively stemming from the two orthogonal contributing terms 
appearing in equation (|T8|). In particular, it is possible to analytically evaluate 
the instant of time, i| a | , at which the system exhibits in the best way such parity 
effect-dependent consequences. Thus, manipulating the atomic velocity in such 
a way that the interaction time between the cavity field and the two-level Ry- 
dberg atom exactly coincides with t\ a \ and measuring the internal atomic state 
immediately after the atom leaves the cavity, the bimodal field is projected into 
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a state having the desired character of a bimodal Schrodinger cat. Indeed 
the state can be written as a superposition of two orthogonal field states, 
\ip) and \x), macroscopically distinguishable in the sense that in the bimodal 
state | f) the field energy is essentially concentrated in the initially empty mode 
whereas in the bimodal state \x) all the energy of the field may be found in the 
other mode. 

2.7.2 Using bimodal cavities for generating single-mode nonclassical 
states 

The interaction between single Rydberg atoms and bimodal high-Q cavities can 
be also exploited with success for generating interesting nonclassical states of a 
single cavity mode. As we will show in this section, the presence of a second 
quantized mode of the electromagnetic field turns out to be very useful playing 
the role of an auxiliary mode whose initial condition may be appropriately cho- 
sen in order to optimize single-mode nonclassical state generation procedures. 

In for example, a bimodal high-Q cavity is used for the controlled prepa- 
ration of single-mode Fock states. 

This method, developed in the context of the nondegenerate two-photon 
micromaser, is based on appropriate sequences of wave packet reductions so 
that its degree of success is accordingly measured in terms of a conditional 
probability. Qualitatively speaking, the main result reached in this article is the 
possibility of conjugating experimentally acceptable values of such a probability 
with a reasonable high value of the number of photons deposited in the cavity 
mode of interest. Following the method presented in Q the generation, for 
example, of the Fock state 1 10) is greater than 0.1. The practical reliability 
of the proposal reported in [p3| is moreover carefully discussed from several 
points of view. In particular the nonideal behaviour of the atomic velocity 
selectors today available, has been taken into account by the authors. Their 
analysis puts clearly into evidence that the procedure is characterized by an 
experimentally significant stability against the unavoidable uncertainty in the 
atom-field interaction times. 

In the context of this topical review, it is important to point out that, as 
shown in [ p3| , the possibility of having experimentally interesting values of the 
probability of success, directly stems both from the peculiar two-photon two- 
mode interaction mechanism adopted in the paper, and, above all, from the 
existence of an auxiliary mode. In the context of their proposal, in fact, the 
initial photon distribution of such a mode is singled out looking for the best com- 
promise between purely mathematical optimization arguments, concerning the 
probability of success of the scheme itself, and practical limitations stemming 
from considerations of an experimental nature. Guided by these considerations, 
the bimodal cavity is prepared leaving the mode of interest in its vacuum state 
and exciting the auxiliary one in a coherent state |a) with intensity \a\ 2 appro- 
priately chosen. Starting from this initial condition and inverting the role of 
the two modes, that is considering as mode of interest the one prepared in the 
coherent state, it is also possible to generate a class of quantum superpositions 
of two single-mode coherent states each one having the same intensity of the 
initial state |9£| In more detail, exploiting the passage of one atom only 
through the bimodal cavity initially prepared in accordance with the request 
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previously discussed, a conditional measurement of the internal atomic state, 
after the atom has left the resonator, projects the mode of interest in a state 
having the following form: 

\*}=A[\p)-e- iXr \/3e- iXT }} (19) 

A being an appropriate normalization constant. In equation ( |l9| ) we have 
indicated respectively by A and r the coupling constant and the interaction 
time between the atom and the bimodal cavity. Moreover /3 = ae" n ^^ and 
7_ (r) = (-| — 0Ji)t, lo\ being the frequency of the mode prepared in the state 
(J19|). In the next section we present a class of dynamical problems studied in 
the trapped ion domain. This class includes some of the problems examined 
in the course of this section reinterpreted in the new experimental context. In 
addition we shall study the dynamical properties of other hamiltonian models 
extraneous to CQED domain, describing the laser driven coupling of internal 
and external ionic degrees of freedom, still in terms of pseudo-spin and two 
bosonic dynamical variables. 
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3 Ions in bidimensional traps 



3.1 Introduction 



Recent progress in laser cooling and trapping techniques has opened up new 
perspectives for research in the fields of atomic physics and quantum optics. By 
cooling down ions to very low temperatures, indeed, purely quantum manifes- 
tations become observable 



100, 101] 



A single ion confined in a Paul microtrap is theoretically equivalent to a 
particle, moving in a harmonic potential, whose center of mass (cm.) motion 
is quantized as harmonic oscillator. By irradiating the ion with classical laser 
beams, its internal (electronic) and motional degrees of freedom can be coupled 
due to the momentum exchange between the laser field and the ion. 

It has been shown |102|, |103|, 104, 106] that there exist suitable conditions 



under which the effective Hamiltonian describing such a system has a form sim- 
ilar to the JCM where the quantized radiation field is replaced by the quantized 
ionic cm. motion. This circumstance suggested the idea of studying phenom- 
ena typical of the cavity QED context with trapped ion systems. A comparison 
of experiments measuring quantum Rabi oscillations in the cavity QED and 
trapped ion contexts [104, |105| shows that the coherence time characterizing 
the dynamics of single trapped ions is at least one order of magnitude bigger 
than that typical of the behaviour of atoms crossing high Q cavities. An impor- 
tant aspect which makes the trapped ion systems so attractive and interesting 
is, indeed, the fact that they are very well isolated from external environment 
being thus very good candidates for experiments relying on the persistence of 
quantum coherences. Moreover, by changing the lasers configuration and pa- 
rameters, a huge class of effective Hamiltonians may be engineered. These facts, 
together with the possibility of measuring the internal ionic state with almost 
unit efficiency, gave rise to many proposals extending Feynman's ideas on the 
use of a given quantum system to simulate the behaviour of another quantum 
system [107|. As far the topic of this paper is concerned, many bimodal inter- 
actions and phenomena typical of both cavity QED and nonlinear optics, can 
be easily obtained in the trapped ion context. For example the vibrational ana- 
logue of parametric amplification, Kerr-type nonlinearities, multimode mixing 
and multiphoton down conversion can be realized [108, 109, |110fl. Moreover a 



Mach-Zcndcr interferometer acting on two modes of oscillation of the trapped 
ion may be simulated allowing the r ealization of entangled states of the ionic 



vibrational 2D motion [141, 112, 113 . Also maximally entangled states and Bell 
states may be engine ered by m eans of several different kinds of easily realizable 
bimodal interactions [112, 113 1. 

It is worth noting, however, that when the spatial wave function representing 
the cm. motion is no longer small compared with the driving laser wavelength, 
there appear nonlinear effects, characteristic of the trapped ion system, which 
have no counterpart in standard nonlinear optics [114, 115]. Under these condi- 
tions interesting phenomena appear due to the overlap of the atomic center of 
mass wave function with the laser waves causing, nonlinear modification of the 
coupling strength gl|, |o| |o|. 

Summing up, trapped ion systems, on one hand allow to verify experimen- 
tally many interesting phenomena which are difficult to observe in other systems, 
and on the other hand they present a very rich quantum dynamics showing non- 
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classical effects having no analogue in optical systems. The considerations we 
have done justify the growing interest toward such systems. In what follows we 
will present a review of the studies on bidimensionally trapped ions. 

We first describe briefly the basic principles of the Paul microtrap and of 
the techniques for detecting the ionic internal states. The effective Hamiltonian 
describing the ionic dynamics is then derived and some proposals for engineer- 
ing generic bimodal interactions are briefly reviewed. In subsection III B we 
concentrate our attention on the dynamics of the laser driven system and, in 
particular, on the appearance of some nonclassical effects in the ionic dynamics. 
Finally, in subsection III C several methods for generating different types of 
nonclassical states of the ion are discussed. 



3.1.1 Quantum motion of ions in Paul traps 

Due to their net charge, atomic ions can be confined by particular arrangements 
of electromagnetic fields. To study ions at low energy, two types of traps are 
usually used: the Penning trap, which uses a combination of static electric 
and magnetic fields, and the Paul or rf trap, which confines ions primarily 
through pondcromotive forces generated by inhomogcncous oscillating fields. 



The operation of these traps is discussed in various reviews [101, 117] and in a 
recent book by Ghosh |l!8| . 

Over the last few years many experimental results for ions confined in Paul 
microtraps have been presented. For this reason we choose to concentrate our 
review to physical situations in which the trapping apparatus is a Paul micro- 
trap. In Fig. ^ we show a schematic diagram of such a trap. If a constant 
voltage U > and an rf alternating voltage V cos(aj r ft) are applied to the trap 
electrodes, an ion of charge q and mass m moving inside the trap sees a time 
dependent potential of the form 

$(z, y, z,t) = q (Q xx (t)x 2 + Q yy (t)f + Q zz (t)z 2 ) (20) 

where x, y and z are the coordinates of the center of mass motion of the ion 
and Qu(t) (i = x,y,z) are the diagonal elements of the traceless tensor of the 
quadrupole moment, given as 



Qu(t) =Q U [U + V cos{u r f t)\ (i = x,y,z) (21) 

From equation ( |20| ) and ( pl| ) the equations of motion of the ionic cm. coor- 
dinates can be derived in the form of Mathieu differential equations. Stable 
solutions of these equations are superpositions of a cm. oscillation at the driv- 
ing frequency u> r f, usually called micromotion, and oscillations of typically much 
smaller secular frequencies Vi in each direction i = x,y, z. Due to the different 
time scales of secular and micro motion, the latter can be neglected by time 
averaging the equations of motion over a period of the r.f. drive. Under these 
conditions, the effective time averaged potential assumes the form: 

$(x,y,z) = ^(vlx* + v 2 y f + v>z 2 ) (22) 

where v x , v y and v z are the secular frequencies in the x, y and z axis of the trap 
and, for quadrupole traps as the Paul trap, satisfy the relation v x +v y — v z . 
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Thus the quantized cm. motion of a single trapped ion in such an oscillating 
electric quadrupole field can be effectively described by a 3D harmonic oscillator 
with the Hamiltonian 



i—x.y,z 



H cm = J2 ^(4^+2) ( 23 ) 



Here a x ^ y ^ and a, x {y,z) are the creation and annihilation operators of vibrational 
quanta in the mode corresponding to the motion along the x{y 1 z) direction, 
respectively. As usual the relation between a*, ^ and a x ( y ^ with the position 
operator x and its conjugate momentum p x is expressible by 



ft /„+ - \ - \rahv x ,„ + „ N 

2^rJ al+ax > P*=*\/—(4-a X ) (24) 

Analogous relations hold for the other two directions y and z. 

3.1.2 Interaction of single trapped ions with external laser fields: 
the JC model 

We first consider the situation where an oscillating classical electromagnetic 
field propagating along the x direction irradiates a single trapped ion which is 
constrained to move in this direction in a harmonic well with frequency v x . In 
more detail we consider a uniform linearly polarized e.m. wave written as 

E(jc, t) = E e 4 ( fcs ^^*) + h.c, (25) 

resonantly driving a single electronic transition between two ionic internal states 
|— ) and |+) . In passing we note that the electric field is an operator acting on 
the Hilbert space of the harmonic oscillator as a consequence of the presence 
of the position operator x in the exponent. The Hamiltonian describing the 
system under scrutiny is the following 

H = H Q + Hint (26) 

with 

Ho = fiv x [a x a x + + hLd a z (27) 

and 

H int = -d • E(i, t) (28) 

where d = d<r + + h.c, d being the complex matrix element of the appropriate 
dipole operator between the states |+) and |— ). 

In the Rotating Wave Approximation (RWA) the interaction Hamiltonian 
assumes the form 

H mt = fiVte l{kS: ' ulLt) a + + h.a. (29) 
where the coupling constant f2 is defined as 



« = -^ (30) 
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Passing to the interaction picture with respect to Hq we get 



It is useful to expand the exponential operator e lkx in its normally ordered form: 



mill 



m-\-l 



..1=0 



In equation ( |32| ) we have introduced the Lamb-Dicke parameter 77 = y 2 m^ • 
This parameter is proportional to the ratio of the center of mass uncertainty 
Ax = y/h/2 mv x in its vibrational ground state and the laser wavelength A = 
2tt jk. The physical situation in which 



+ 1 < 1, (33) 

where (n x ) is the mean number of oscillatory quanta along the x direction, is 
commonly referred to as the Lamb-Dicke limit. In such a condition the ion 
couples with the laser behaving as a point-like object. Alternatively, rj can be 
seen as the ratio between the photon momentum fik and the ionic ground state 
momentum uncertainty Ap x = yJfimv x /2. Thus, in the Lamb-Dicke regime, 
the effect on the center of mass wavepacket due to momentum transfer during 
absorption or emission of a photon is negligible. We wish to underline that very 
often in literature the condition 77 1 is referred as Lamb-Dicke limit. This 
condition is less restrictive than the one expressed by inequality (|33|). In the 
following we will also refer to r\ -C 1 as the Lamb-Dicke condition but having in 
mind that the mean number of oscillatory quanta have to be small enough to 
satisfy ©. 

Inserting equation (p2|) into equation(|3l|) yields 



H( nt (t) = hn<j + e-v / 2 £ i^^at^e^-^™-^]' + h.c. (34) 

m,l=Q 

This Hamiltonian consists of a sum of terms oscillating at frequencies dif- 
fering from loo — lul by multiples of the vibrational frequency v x . In the so 
called resolved-sideband regime defined by the condition v x 3> O, L, P being the 
total linewidth of the electronic transition, the different vibrational sidebands 
are clearly resolved and it is possible to tune the laser in resonance with a spe- 
cific vibrational sideband. Moreover, in such a physical situation, if the laser is 
tuned to the fcth sideband, with k > and k < corresponding to red and blue 
detuning respectively, 

u L = ujq - kv x , (35) 

we may neglect in equation (|34|) all terms oscillating with multiple vibrational 
frequencies. As a consequence the effective interaction hamiltonian in the inter- 
action picture looses its time dependence, becoming 

Hf nt (t)~H? nt = nn& + f k (ata x ,r,)(i V a x ) k + h.c. k>0 (36) 

HL(t) =s H&t = (w?4r fe f-h (4m) + h - c - k < 
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The operator valued function fk {a\.a x ,r]) , with k > and in its normal ordered 
form, can be written down as follows [114 



1=0 



{l + k)H 



allaL 



k > 



(37) 



The interaction Hamiltonian given by equation (|36|), due to the presence of 
the operator ff. (a x a x ,'qj, represents a nonlinear k-quantum Jaynes-Cummings 
model. There exist conditions under which this Hamiltonian reduces to the 
well known JC Hamiltonian and multiquantum JC Hamiltonian. To better 
understand this point consider the vibronic Rabi frequency defined as 



Q n ,n+k = fi(n|/k (al& x , rj) (irja x ) h \n + k} = (ir)) k tte 1,2/2 1 



(n + k) 



■4 fc) (v 2 ) ■ 

(38) 

Equation ([58|) shows a nonlinear dependence on n departing from that char- 
acterizing the usual /c-quantum JC model. Moreover, due to the presence of 
the associated Laguerre polinomials in equation (|8|) , there exist pairs of values 
of r\ and n > in correspondence of which the laser-atom coupling practically 
vanishes. Let us consider the nonlinear k — 1-quantum JC interaction. In 
the Lamb-Dicke regime, the operator f% appearing in the Hamiltonian can be 
simplified as f\ ~ 1 and thus the vibronic Rabi frequency assumes the form 
f2 ni „_l_i ~ irj$l\fn + 1 showing the typical behaviour of the standard JC model. 
The same considerations hold for a generic value of k leading us to the follow- 
ing conclusions. The light mediated interaction between the quantized motion 
of the cm. of an ion moving in a Paul microtrap and its electronic degrees 
of freedom shows a nonlinear character. This behaviour, traceable back to the 
specific adopted vibronic-coupling mechanism, is typi cal of thi s system and does 
not have analogues in cavity QED quantum models [114, 119 . However, in the 



Lamb-Dicke regime, that is when the cm. wavepacket is well localized in com- 
parison with the laser wavelenght, the nonlinear multiquantum JC interactions 
reduce to the multiquantum JC models common in quantum optics. Nonethe- 
less, it is worth stressing that ionic nonlinearities play in general, an important 
role in the experiments and very often need to be taken into account to describe 
correctly the dynamics of the system[104, |112fl . 

We wish to underline that many experiments have been performed, during 
the last few years, verifying the theoretical description of the ID Hamiltonian 
model given by equation^) [104]. Moreover, a huge class of vibrational states 
such as Fo ck st ates, coherent states, squeezed states [104 and Schrodinger cat 
like states [120 , have been experimentally realized and even the tomographic re- 
construction of the Wigner function of a vibrational Fock state has been demon- 
strated 1 121]. In most of the experiments the electronic states of the ion are 
coupled by using two photon stimu lated Raman t rans itions through a third 
optical level as indicated in Fig. [?] |104 |120| , |I12[ |12]| . For sufficiently large 
detunings A such that the third intermediate state can be adiabatically elim- 
inated, only two electronic levels are involved in the dynamics. In this sense 
stimulated Raman transitions are formally equivalent to a narrow single photon 
transition 104 , |112| . Moreover they recombine the advantages of strong opti- 
cal electric field gradients (allowing manipulation of the state of motion) and 



2.x 



microwave stability of the crucial difference frequency. It has been shown that 



[104, 112 1 the effective Hamiltonian model in this case has the same form of 
the one given by equation (p6[), where now the coupling constant is given by 

d .E <1<2>) 

£1 = Qif^/A, with r^i(2) = 1<2) h and A detuning from the third optical 

level, and i] = \j k = ki k2. We indicate with di and d.2, 

(2) 

and Eq , ki and k2 the dipole matrix elements, laser amplitudes and wave 
vectors of the two lasers guiding the Raman transition respectively. As a con- 
sequence, the Lamb-Dicke parameter, depends, in this case, also on the relative 
angle between the two incident lasers through the modulus of k. In this way it 
is possible at the same time to excite the ion motion along one direction only 
and to manipulate easily the L-D parameter simply by changing the relative 
directions of the beams irradiating the ion. 

Raman type configurations can be employed also to implement vibrational 
transitions without affecting the electronic state of the ion as shown in Fig.^j. 
The form of the effective interaction Hamiltonian in the interaction picture is, in 
this case, very similar to the one given by equations (|3^) and (^), but electronic 
inversion operators are absent: 

Hf n t = mf k (ata x ,r))(ir)a x ) k + h.a. (k > 0) (39) 

Htt = hSl(ir]aiy k f- k (alax,v)+h.a. (fc < 0) 

We conclude this subsection sketching briefly a very efficient method, based 
on the so called quantum jumps technique, used in the experiments for measur- 
ing the population of the electronic states of the ion ]122|] , It is worth noting 
that such a quantity is the only one directly measurable for trapped ion sys- 
tems. The method essentially consists in coupling the ground electronic level 
|— ) to a third auxiliary level |r) by means of a third laser beam, as shown in 
Fig|| The transition |— ) |r) is chosen to be a dipole allowed one. In this 
conditions, the presence of fluorescence detects the ion in the electronic ground 
level and its absence in the excited level. In other words, the detection of null 
fluorescence prepares the ion in the excited state since after the measurement 
the quantum system collapses in the eigenstate correspondent to the measured 
outcome. With this method also if a modest number of scattered photons are 
detected, the efficiency of the ability to discriminate between the two states 
approaches 100%. 



3.1.3 Generic bidimensional Hamiltonian models 

In this subsection we wish to present the extension of the ID Hamiltonian model 
given by equation ( |36| ) to the case wherein two independent modes of motion 
of the ionic cm., in the effective harmonic potential present inside the trap, are 
excited. In the previous subsection we have assumed that the laser interacting 
with the ion propagates along the x-axis of the trap. When the laser wave vector 
lies on a plane of the trap, say the x-y plane, the dynamics couples the motion 
on that plane only. In this context, usually referred to as 2D trapped ion, a 
derivation analogue to the one given in |l08| leads to the following expression 
for the general effective Hamiltonian in the interaction picture 

H ZV mv = h^e^c/mx (4>«*,»7x) <?m y (a\,a y ,r) y ) a\ + h.c. (40) 
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with 



(«j>%:%) = (iVjVj) f\k\ (ajaj,Vj) > 0) (41) 

('V '';•'/..) = f\k\ (a]a>j,Vj) (iVjaj) lkl (k < 0) (42) 



,9fc 

3fc 



where we have assumed the resonance condition cj^ — cl>o = tti x v x + m y v y and we 
have used the property = v x +v y , valid for all quadrupole trapping potentials. 



In equations (41) and ( |12| ) the operator function is defined by equation (p7|). 
The parameter e is equal to 1 when internal transitions are involved (Fig.Q) and 
e = when the internal state is unchanged (Fig. ||). In the Lamb-Dicke regime, 
that is assuming r\ x ,r\ y << 1, one gets 

H™f°> m « <0 = hn' (a*) 1 ™* 1 {a y f^ o% + h.c. (43) 



L int — V*VJ " + 

J?™f 0,m ' <0 = ftfi' (4)' m *' {aj m ^ a% + h.c. (44) 



#mi = K)' K) cr + + /l.C. (45) 

^> , mH >o = m , (4 y-i (a ty«.i ^ + h c (46) 

(47) 

with f2' = n(ir?a;)' mx ' (z77j,)' my ' . Hamiltonian models similar to the ones ex- 



pressed by equations (j4J) and (45) are considered in |10q , 123 1. Both in [123 1 
and in [ 108 1 it is shown that, in the case e = 0, an appropriate choice of the 
detunings of the lasers leads to effective Hamiltonians typical of the non linear 
optics context. We will analyze in more detail some examples in Section II. 

It is important to underline that the derivation of the effective Hamiltoni- 
ans ( p0[ ) with the methods presented in 1 108 , 123| , 112, 124 is valid under the 



assumption that the trap frequencies v x , v y and v z are incommensurate. When 
this condition is violated, as for example in the case of bidimensional isotropic 
traps for which v x = v y , "unwanted" resonant terms may appear in the effective 
interaction Hamiltonian in the interaction picture. For isotropic traps, in order 
to avoid the appearance of such terms, it is necessary to use more complex laser 
configurations. 

3.2 Dynamics 

3.2.1 Trapped ions as quantum simulators 

In this section we describe the dynamics of some 2D trapped ion models studied 
in literature over the last few years. As we have already pointed out, initial 
studies on single trapped ions put into evidence the formal analogy existing, 
under certain conditions, between these systems and some nonlinear optical 
systems. As a consequence, when these physical conditions are satisfied, the 
Hamiltonian models given by equations ([l3|)-([l6]), can lead to simulations of 
various devices of practical interest. As an example, a (50/50) beamsplitter can 
be simulated. The action of the beamsplitter on vibrational Fock states has 
been studied by Gou and Knight in the case of 2D isotropically confined ions 



[ 1 1 1 1 . The authors show that for these systems the unitary operator describing 



the beamsplitter coincides with that one describing rotations of the x — y axes 
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of the trap around the z axis in the real space. Similarly in [113, 112| a Mach- 
Zender interferometer acting on the two modes of oscillation of a single trapped 
ion, has been considered. The analogy with the Mach-Zender interferometer for 
bosons is that the two input modes to the boson interferometer are replaced by 
the x and y modes of ion oscillation. 

Another example of the use of trapped ions to simulate the dynamics of 



optical sistems is given in [110]. In this paper, indeed, Agarwal and Banerji 
demonstrate that a single bidimcnsionally confined ion provides a natural way 
of realizing the quantum limit of second harmonic generation in a degenerate 
optical parametric oscillator (OPO). They show that simulating the quantum 
regime of OPO allows to demonstrate the existence of an effect strictly quantum 
in nature. 

3.2.2 Motional nonlinearities 

All the examples concerning the use of trapped ions as quantum simulators rely 
on the assumption that the Lamb-Dicke condition is satisfied. However, when 
the spatial extension of the atomic wave function representing the cm. motion 
is no longer small compared with the driving laser wavelenght, nonlinear effects 



having no counterpart in standard nonlinear optics emerge. In [105, 119| several 
Raman-type excitations for inducing various kinds of nonlinear interactions in 
the quantized motion of a trapped ion are presented. It is there shown that stan- 
dard effects can be realized, including coherent displacement, Kerr nonlinearities 
and parametric modes couplings. In the laser-assisted motional dynamics, how- 
ever, far from the Lamb-Dicke regime, additional nonlinearities appear because 
of the effects on the atomic wave function due to the driving laser waves. One 
of the consequences of these nonlinearities consists in a nonlinear partitioning of 
the phase space, the action of the time evolution being different in neighboring 



phace-space zones [108]. To better understand the origin of the phase space 
partioning let us consider first, for semplicity, the single mode dynamics, where 
only the motion in the x direction is affected by the lasers (rj y = rj z = 0). In 



this case the Hamiltonian is given by equation (f39|) with fk (a x a x ,r]) given by 



equation (37). If the Raman lasers are tuned to the first vibrational sideband, 
ujl = ijJli — UJ L2 = Vx, the structure of the time evolution operator shows some 
formal resemblance to a nonlinearly modified coherent displacement operator: 

Uint(t) = D -rj x Cl*fi(h x ,ri x )t = exp -rj x O.*talfi(n x ,rj x ) + h.c. (48) 

Note that for small values of the Lamb-Dicke parameter, r) x -C 1, the operator 
( p8[ ) may be replaced by the usual displacement operator D(—r] x Vl*t/2). Insight 
into the nonlinear modifications of the displacement can be obtained by con- 
sidering the expectation value of the operator function fi(n x ,f] x ) in a coherent 
state | a): 

(a\f\(h x ,r] x )\a) = (vM)^ -h^M)^ 1 ' 2 (49) 

This function shows the oscillatory behaviour typical of the Bessel function, 
where the zeros are of special importance. Indeed these zeroes lead to the 
partitioning of phase-space into different (radial) zones, with their boundaries 
being defined by the zeros. In neighbouring phase-space zones the actions of the 
time evolution appear to be significantly different from each other. For example, 
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one can see from numerical calculations that in two adjacent zones a nonlinear ly 
modified "displacement" operator acts in opposite directions. Consequently, a 
quantum state initially located on the boundary between two zones will be split 
in two subspaces which eventually give rise to quantum interferences. 

When two or three dimensions are involved in the Raman induced dynamics, 
the partioning effect appears in the phase space of each motional degree of free- 
dom. As a consequence, the coupling between different motional modes will be 
strongly influenced by the inter play of these nonlinear effects. So, as pointed out 
by Wallentowitz and Vogel in flO§[| , besides the feasibility of realizing phenom- 
ena well known from linear optics in the motion of bidimensionally trapped ions, 
these systems open up novel possibilities for studying new nonlinear quantum 
phenomena. 



3.2.3 Quantum effects in the dynamics of bidimensionally confined 
ions 

In this subsection we discuss some quantum effects recently brought to light in 



the dynamics of 2D isotropically confined ions [126, 127]. Quantum phenomena 
similar to the ones we are going to describe have been also studied in cavity 
QED context. Nonetheless, in the context of trapped ions they acquire a new 
physical meaning and they appear to be much easier to verify experimentally. 
Before describing the system under scrutiny we want to underline that here, as 
well as in many of the theoretical proposals we will describe in this review paper, 
an optical transition coupling directly the electronic states is considered. In the 
experiments, however, two photon Raman transitions are preferred. However, it 
is always possible to substitute the direct coupling scheme with the two photon 



Raman scheme obtaining the same effective Hamiltonian models [ 104 1 . 

We consider a two-level ion confined in a bidimensional isotropic harmonic 
potential characterised by the trap frequency v x = v y = v. It is well known 
that the bidimensional harmonic oscillator may be associated to a generalised 
Schwinger angular momentum operator J = (Ji, J2, J3) as follows 

t n.. -4- r7.t/7._ . nJ a... — n ~ . /7.t/7._ — nJri 



a x a y + a y a x ~ _ a x a y - a y a x - _ a x a x 



Jl — ^ J2 — TT. J3 — 



ClyCly 



2 2i 2 

(50) 

It is worth stressing that these generalized angular momentum operators ac- 
quire, in the context of 2D isotropically confined ions, a direct and interesting 
physical meaning. The operator J3 is indeed simply half the difference between 
the number of vibrational quanta along the x and y axis of the trap. The 
operator J2 is proportional to the z-component of the ionic adimensional angu- 
lar momentum operator L z = j- (xp y — yp x ) = —i (a x a y — a y a x ) — 2^2- Finally 
the operator J\ has the property of being proportional, under certain conditions, 
both to the quantum covariance C{x,y) — (xy) — (x}(y) between the position 
operators x, y and to the quantum cova rianc e C{p x ,p y ) — (p x Py) — (Px)(Py) 
between the momentum operators p x , p v [|l26[ . 

It has been shown that by irradiating the ion with an appropriate configu- 
ration of laser beams the physical system can be studied, in the L-D limit and 



in the interaction picture, by the following Hamiltonian model [144, |126|, 127 



129, 128, 130 



H int = Ml' (a x a y a + + h.a.) (51) 



32 



with f2' = fir] 2 . This Hamiltonian is unitarily equivalent to the last term of 
equation ( |l6| ) responsible for the parity effect presented in section II D 4, in 
the context of cavity QED. In this case, however, Stark shift terms present in 
equation are absent and the parity dependent quantum effect discussed 
in section II D 4 is much "cleaner". Moreover it turns out that new physical 
observable, having a direct physical meaning only in the trapped ion context, 
show an interesting nonclassical behaviour. In particular in ]I26| it is shown 
that, starting from an initial highly excited unidimensional Fock state \n s — 
N,riy — 0) with N 3> 1, along the direction x having an angle of 7r/4 relative 
to the x axis, the time evolution of the mean value of the operator J\ may be 
written as follows: 

<J l( i)>^ CO sW^) (52) 



2 V N , 

From equation (B2) it is easy to deduce that at the time instant t/v = tjt one 



n* 

gets (Ji(t)) — (— I)^ -1 -^- This result clearly evidences that the system under 
scrutiny possesses an inherent peculiar sensibility to the parity of the initial total 
number N of vibrational quanta. There exists indeed a A-dependent instant 
of time ijv at which (Ji(t)) assumes values crucially related to the evenness or 
oddncss of N. In more detail, ( Ji(tiv)) is equal to —N/2 when N is even and 
+N/2 when A is odd. Fig.[l(] shows the time evolution of the quantity (Ji(t)). 
Remembering the link between (Ji(t)) and the quantum covariances C(X,Y) 
and C(P X ,P V ) and taking into consideration that the maximum (minimum) 
mean value of J\ is N/2 (—N/2), this result can be readily interpreted in the 
following way. The dynamics of the system under scrutiny is characterised by the 
existence of a N-dependent instant of time at which the ionic oscillatory motions 
along the x and y axis are maximally or minimally correlated in position and 
momentum in dependence on the oddness or evenness of N respectively. 

Another interesting quantity displaying a parity dependent temporal evolu- 
tion is the variance of the operator J2 which is proportional to the variance of 
the axial orbital angular momentum operator L z . In a system identical to the 
one previously described, the mean value of J2 is null at any time instant t, 
while at the time instant tjv /2 = nN/2fl the square of the variance assumes the 
value fL27| 

N 2 

(AJ 2 ) 2 (t N /2) ~ (l - (-1)*) — + 0(N) (53) 

This means that (A^) 2 ^^ /2) is of the order of TV ^ 1 when the initial number 
of vibrational quanta is even and of the order of A 2 when it's odd. In other 
words the variation of one vibrational quantum only in the initial conditions of 
the system, quite drastically modifies the dispersion around zero of the measured 
values of the ionic angular momentum z-component at a subsequent specific 



instant of time. In [129 1 we have shown that at this time instant, tjv/2, the 
vibrational and electronic degrees of freedom of the ion are disentangled for odd 
N and maximally entangled for even N. At t = fjv/2, indeed, the probability 
of finding the ion in its ground electronic state is 1/2 or 1 in correspondence 
to N even or odd respectively. We have also checked the robustness of this 
parity dependent entanglement effect against the non-dissipative decoherence 
due to laser intensity fluctuations and imperfect initial state preparation finding 
that, for values of the laser parameters currently used in the experiments, the 
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nonclassical effect, althought attenuated, is still present [131 



132] 



3.3 Generation of nonclassical states 

3.3.1 Synthesis of arbitrary states of two-dimensional vibrational 
motion 

As we have already mentioned in section III A, trapped ions are exceptionally 
well suited systems to study fundamental aspects of quantum theory. By ir- 
radiating the ion with laser fields coupling its internal and external degrees of 
freedom, a single ion can be cooled to the vibrational ground state of the trap 
and its motion can be manipulated coherently to generate nonclassical states of 
motion. The ability to synthesize arbitrary motional states and to control their 
evolution not only allow to verify experimentally fundamental concepts of quan- 
tum mechanics, but also opens new horizons for potential technical applications 
such as quantum information processing. For these reasons, during the last few 
years, several proposals for engineering arbitrary nonclassical states of the ID, 
2D or 3D oscillatory motion of the ionic center of mass have been presented (see 



[119 1 and references therein for a review). In what follows we concentrate our 
attention on bidimensional engineering schemes. 

The first scheme for preparing an arbitrary motional state of a bidimen- 



sionally trapped ion was proposed by Gardiner, Cirac and Zoller [133] as a 
key prerequisite for quantum measurement of an arbitrary motional observable. 
The method, conceptually similar to the one proposed by Law and Eberly in 



the context of cavity QED [ 134 1 , consists in preparing a two- mode target state 
of the form 

i*w> = E E ), (54) 

n x =0 riy—0 

starting from the initial vibrational and electronic ground state \n x = 0,n y — 
0, — ). Such a state is generated by means of an appropriate sequence of laser 
pulses used to coherently distribute the initial amplitude to form a desired 
arbitrary superposition. 



Although conceptually elegant, the method proposed in [133] is difficult to 
implement since the number of laser pulses necessary for preparing the bimodal 
target state depends exponentially on the dimensionality of the subspace of the 
Fock space in which the target state is embedded. Ideal synthesis of motional 
quantum states assumes indeed negligible dissipation effects and perfect control 
of areas and phases of the laser pulses and switching times. Fluctuations of these 
quantities cause some technical noise affecting the efficiency of the schemes. For 
these reasons it is desirable to obtain preparation schemes based on a number 
of laser pulses as small as possible. In this spirit, during the last few years, 
novel approaches allowing to reduce consistently the number of operations re- 



quired for synthesizing desired target states have been proposed [135, 137, 138 1. 
In particular very recently Zheng presented a method for preparing arbitrary 
bidimensional vibrational states of a trapped ion of the form (pj) requiring only 



(M m ax + 2)(N max + 1) laser pulses [138]. It is worth noting, however, that a 



state of the form ( p4[ ) involves (M max + l)(N max + 1) complex coefficients, each 
controlled by the amplitude and phase of a laser pulse. Thus it should exist, in 
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principle, a better approach requiring only (M max + l)(N max + 1) operations 

All the procedures we have described till now concern the generation of ar- 
bitrary bimodal oscillatory states of the cm. motion of an ion confined in a 
trapping potential. They rely on the possibility of engineering suitable interac- 
tion Hamiltonians, generally involving the two mode bosonic operators and the 
atomic ones, to control the evolution from the initial state to the desired target 
state. It has been shown that, like in the case of bimodal cavity fields, such 
bimodal interactions can also be used to synthesize more efficiently arbitrary 
unidimcnsional vibrational states. In a recent proposal, Hladky, Drobny and 
Buzek show how to engineer an arbitrary unitary operator V , transforming any 
state |$) of a single bosonic mode, e.g. a ID vibrat iona l state of a trapped ion 
along the x direction, to another state |^') = V\^} |]l40| . The realization of the 
unitary transformation V is based on an enlargement of the Hilbert subspace 
of the system. In the case under scrutiny the ionic cm. mass mode of oscil- 
lation along x becomes entangled with auxiliary degrees of freedom (ancilla) 
which are represented by the vibrational mode along the y direction and three 
internal electronic levels of the ion. The physical realization of the desired op- 
erator V consists again of a sequential switching of laser fields irradiating the 
ion, the correspondent interaction Hamiltonians being typical and quite easy to 
realize in the trapped ion context. The overall algorithm involves a conditional 
measurement of an internal ionic state, thus the process is not unitary. It is, 
however, universal since the probability of the conditional measurement does 
not depend on the initial vibrational state. 



3.3.2 Bimodal generalized coherent states 

In this subsection we will review some recent proposals for realizing, in the 



context of trapped ions, generalized coherent states (GCS)[141| such as pair 
coherent states, SU(1,1) coherent states and SU(2) coherent states. Over the 
last several years there has been much discussion on GCS of the quantized 
electromagnetic field since they exhibit strong nonclassical properties such as 
sub-Poissonian statistics, squeezing, violation of Cauchy-Schwarz inequality and 



so on 1 145 . Unfortunately they are not very easy to realize in practice in the 
optical context. We have already seen, however, that it is possible to imagine 
rather simple mechanisms by which a great variety of effective Hamiltonians can 
be engineered in the trapped ion context. Such interaction models can be used, 
as we shall see in the following, to generate several classes of vibrational states 
including generalized coherent states. There exist basically two main ways to 
generate a given vibrational state with trapped ions. The first way consists in 
the coherent evolution due to an appropriately engineered interaction Hamilto- 
nian, starting from a given initial state. Eventually a conditional measurement 
can be performed to project the ionic vibrational state in the desired target 
state. The second way takes into account also the effects of decoherence due to 
spontaneous emission in the process of state generation. An advantage of this 
last method is that the desired states are obtained as steady state solutions of 
the master equation containing dissipative terms due to spontaneous emission. 
These target states are called dark states because their generation experimen- 
tally corresponds to the extinction of fluorescence by the ion. Whenever a given 
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class of states is denned by an eigenvalue equation, it is possible to imagine, 
in the trapped ion context, a configuration able to realize the dark target state 
starting from a proper initial condition of the system. To better understand 
this point we remind that, as the damping of vibrational quanta can be signifi- 
cantly suppressed in an ion trap, the dominant decay process is the spontaneous 
emission from the two-level ion. Thus the time evolution of the system, in the 
interaction picture, can be described by a density operator p obeying the master 
equation 



dp 
~dt 



h! 



int ) P 



(2cr_p<T + — o+o_p — pa+a. 



(55) 



where T is the spontaneous decay rate of the excited ionic state and g is the 
modified density operator 



f 



du dv W(u, v y k ^z+vy) p e ^Hux+vy) 



(56) 



accounting for the momentum transfer in the x — y plane due to spontaneous 
emission. In equation ( |56| ) W(u, v) describes the angular distribution of the 
spontaneous emission. The steady state solution of equation (j55|) has the form 



ps = where ps satisfies the restriction 

interaction Hamiltonian has the form 



H Lt>Ps 



0. If the 



.4 



vib 



cr + + h.a. 



(57) 



, where A v a is a generic vibrational operator, then the stationary state of the 
system will satisfy the eigenvalue equation H( nt \^) = el*). 



a) SU(2) coherent states 

We have seen in subsection III B 1 that a Mach-Zender interferometer can be 
simulated with 2D trapped ion systems. An SU(2) interferometer is a peculiar 
kind of Mach-Zender interferometer of particular relevance in quantum optical 
theory. The vibrational SU(2) states are anticorrelated bimodal states formally 
identical to the photon states on the output ports of a SU(2) interferometer with 
number states inputs. As a consequence a way to generate vibrational SU(2) 
states is to prepare an initial bimodal Fock state and then to realize an effective 
Hamiltonian "simulating" the action of a SU(2) interferometer. In Jl 1 1| it is 
shown that the mode-mixing process on the output ports of the interferometer 
is replaced, in the context of trapped ions, by a real-space rotation in the x-y 
plane. This means, in other words, that a bimodal Fock state excited along the 
x — y axis is an SU(2) state in the coordinate system rotated of an angle 9 with 
respect to x — y. We wish to remind that, as seen in subsection III B 3, an 
isotropic bidimensional harmonic oscillator may be associated to a generalized 
Schwinger angular momentum operator J whose components are defined in 
equation (50). The two mode Fock state \n x ,n y ) can be described as a pseudo 
angular momentum state \j, m) which is the common eigenstate of the angular 
momentum operators J 2 and J3 with correspondent eigenvalues j = \ {n x +n y ) 
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and m = \{n x — n y ). The generalized SU(2) coherent states are defined by the 
action of the generalized SU(2) displacement operator as follows 

\T,j)=e(f >S +-P*-)\j,-j) (58) 

where J+ = J x + iJ 2 , J- = Ji - Ui, P = ae^^ (0 < a < vr, < 7 < 2vr) 
and r = tan(a/2)e' _l7 - ) . Reminding that the operator J2 is proportional to the 
z component of the orbital angular momentum operator generating rotations in 
the x — y plane, it is not difficult to convince oneself that the method proposed 
in [111 I for generating SU(2) states allows also to engineer "real"SU(2) coherent 
states, that is states \r,j) with r real. 

Very recently another method for generating dark SU(2) states of the bidi- 
mensional oscillatory motion of a single trapped ion has been proposed [ 143 1 . 
Such states are eigenstates of both the total number of vibrational quanta and 
the vibrational operator A v n, = a x a y + h.a.. Due to the degeneration of their 
eigenvalue spectrum, however, the dark states obtained as steady state solu- 
tions of the master equation (|55|) would be incoherent superpositions of SU(2) 
states. It is shown in [143] that, for large values of the L-D parameter, motional 
nonlinearities come into play removing the eigenvalues degeneration. In these 
circumstances it is possible to show numerically that the steady state solutions 
approximate well SU(2) states [ 143 1 . 



b) Dark pair and SU(1,1) coherent states 

Let us now discuss two similar methods for generating dark pair coherent 



states and SU(1,1) coherent states 144, 145 1. We firstly review the procedure 
for engineering pair coherent states 144 



The pair coherent states |£; q) are defined as the eigenstates of both the pair 
annihilation operator a x a y and the number difference Q = n x — h y and, in the 
Fock state basis have the form: 

00 A 

lfw)=iV,g- 7|r |=| i + < ! ,<) (59) 

where N q — [|£l~' ? ^g(2|£l)] 1S the normalization constant (I q is the modified 
Bessel function of the first kind of order q). It is then clear that, by engineering 
an Hamiltonian model of the form 

H int = HQ [a x a y — £] <x + + h.a. (60) 

the dark state obtained is an eigenstate of a x a y with eigenvalue £. This means 
that, since the number difference operator Q is a constant of motion in processes 
involving simultaneous pair annihilation or creation, starting from an eigenstate 
of Q, like |<J/o) = |<7, 0), one obtains as dark steady state the pair coherent state 



A similar technique is used in [ 145 1 to create one mode and two mode SU(1,1) 
intelligent states of the vibrational ionic motion. SU(1,1) intelligent states are 
a class of SU(1,1) coherent states defined by the property that an SU(1,1) un- 



certainty relation is equalized [147]. It can be shown |148j that these states are 
solutions of the following eigenvalue problem 

(aa x a y +j3aial)\$)=(\$) (61) 
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Moreover, they are eigenstates of the number difference operator Q. Thus, 
preparing the system in an eigenstate of Q and irradiating the ion with a con- 
figuration of lasers realizing the Hamiltonian model 

HL = h [m x a v + (laial - (] a+ + h.a. (62) 

the steady state solution, indicated by extintion of fluorescence emission by the 
ion, is an intelligent SU(1,1) state. 



c ) Generation of pair coherent states via QND measurements 

To conclude this subsection we wish to describe a recent alternative tech- 



nique for generating pair coherent states [149, 150 1 . The scheme is based on 



a Quantum Non Demolition (QND) measurement procedure and allows actu- 



ally to engineer several classes of vibrational states 151 . The ion is irradiated 
by two laser beams propagating along the x and y axis of the trap and realiz- 
i ng, i n the L-D lim it and in the interaction picture an interaction of the form 



[|4g [150[ [15f| |l52| 



H Lt = l^Lx - &L y - xtyx - n y )\a + + h.a. = \Cl Lx - Q. Ly - xQ\°+ + h.a. (63) 

where Qlx = dE x /h, fl^y = dE y /h and we have assumed the condition Qlx'Hx = 
riLyT]y = x, with Hlx &Ly The nonlinear JCM of equation ( |63| ) is exactly 
solvable and it is not difficult to prove that, if the system is initially prepared in 
a state of the form \n x , n y , — ), then Rabi oscillations of the electronic occupation 
occur, the vibrational state remaining unchanged: 

= cos [Vt(n x ,n y )t] \n x ,n y , -) - isin [Q(n x ,n y )t] \n x ,n y , +), (64) 

with £l(n x , n y ) = — &>Ly ~ X<j\ and q = n x — n y . We have seen in section 
III A that, in a typical trapped ion experiment, a measurement of the internal 
electronic state is realized coupling the ground state with a third level by an 
auxiliary laser beam |I22fl . If a null fluorescence is detected then we may claim 
with certainty that the ion is not in its ground state. In other words, this 
detection process prepares the ion in the excited state. Suppose now that at 
t = the ion is prepared in its electronic ground state whereas the vibrational 
motion is in a generic superposition of two-mode Fock states. Equation ( p4[ ) 
suggests a simple protocol to generate a prefixed bimodal Fock state \m x ,m y ) 
provided that the probability of finding at t = the center of mass motion in 
the state \m x ,m y ) is different from zero. In order to generate this state, we 
imagine to measure several times the internal electronic state at time instants 
ti appropriately chosen. 

The first measurement is performed at the time t\ which satisfies the con- 
dition £l(m x ,m y )ti = (2Zi + l)vr/2, (h — 0,1,2,- ■ •). If the null fluorescence 
is detected, according to equation (|64|), the vibrational state \m x m y ) keeps 
the original probability unchanged whereas the other components present at 
t = decrease their probabilities. In other words, after the renormalization of 
the state, the probability of the state \m x m v ) is increased with respect to the 
other two-mode Fock states. This process can be obviously iterated measuring 
the null fluorescence at time instants tk (k = 2, 3, • • •) defined by the condition 
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Q(m x , m y )tk — ifc7r, (Ik = 1, 2, 3, • • •). It is easy to convince oneself that, follow- 
ing this procedure, the prefixed state|m x ,mj,) can be generated. Remembering 
the Fock state expansion of the pair coherent states, given by equation (|59|), 
and noting that all the states |n + q, n) share, in the system under scrutiny, the 
same Rabi frequency Ct(m x — n + q, m y — n), it is possible to persuade oneself 
that the procedure described before may be used in our context to engineer pair 
coherent states. It turns out that, starting from a bimodal coherent state \a, (3) 
it is possible to generate the pair coherent state |£ = a(3,q) with fairly good 
efficiency. As an example, Fig.[y] shows the probability distribution relative to 
the vibrational state \a/3,2) in correspondence to \a\ 2 = 2.5 and \f3\ 2 = 1.5, 
and x/{^Lx — ^Ly) = 0.0049. It is worth stressing that this quantum state 
manipulation procedure requires several successful null fluorescence detections 
in order to generate the desired state. In the case of the example of Figjll] the 
probability of success is P = 0.172, which is already a significant value from an 
experimental point of view. 



3.3.3 Schrodinger cat-like states 



Since the famous cat thought experiment suggested by Schrodinger in 1935 [ 153 j , 
superpositions of macroscopically distinguishable quantum states, also known as 
Schrodinger cat states, have become a longstanding example of the peculiarities 
which occur in the interpretation of quantum reality. To explore and gain insight 
into the fundamental issues of quantum theory, a number of schemes have been 
proposed for the realization of such states. However, the generation of cat states 
and the detection of their quantum properties is an extremely hard tool since 
decoherence effects due to the interaction with the external environment tend 
to transform these quantum superpositions into statistical mixtures destroying 
their coherence properties. For this reason single trapped ions have shown to 
be very well suited to the generation of Schrodinger cat states, since, as we have 
already pointed out, dissipative effects which are inevitable from cavity damping 
in the optical or microwave regime, can be significantly suppressed for the ion 
motion due to the extremely weak coupling between the vibrational modes and 
the external environment. The archetype of a Schrodinger cat states is given by 
the superposition 

\9)cat = N(\a) + e*\-a)), (65) 

where a is a coherent state of a single mode quantized field and N is a nor- 
malization constant. These states are referred to as even, odd and Yurke-Stoler 



coherent states when (ft — 0, ir and 7r/2 respectively [ 154 1 . Schrodi nger cat states 
of this type have been recently realized for a trapped y £?e + ion |12C|. With an 
average of about nine vibrational quanta, wave packets of maximum separation 
of about 83 nm, significantly larger than the size of a single component wave 
packet, about 7 nm, were obtained. 



a) Bimodal Schrodinger cat states 

In this subsection we review some proposals for generating Schrodinger cat- 
like states of the bidimensional oscillatory motion of a single trapped ion. The 
first generation scheme we consider 1 152 allows to engineer cat-like states which 
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are the generalization to two modes of motion of the states experimentally real- 



ized in [120]. The system considered is an isotropically confined ion irradiated 
by laser beams in a configuration similar to the one described in the previous 
paragraph and realizing the effective Hamiltonian ([33]). The only difference 
is that here the Rabi frequencies of the two lasers are assumed to be equal 
Qlx = &Ly = With this restriction, the interaction Hamiltonian in the 
interaction picture assumes the form 

Hl nt = ~x(n x - fiy) (o+ + ctI) (66) 



As shown in [[152|, if the initial state is a bimodal coherent states \a) = \a, (3) 
and the ion is in its ground internal state, after a time t > the state of the 
system evolves to 

|tt(t)) = |-)|S+) + |+>|S_) (67) 



where 

15+) = 



i (\ae^ 2 ,(3e-^ 2 ) ± jae"^/ 2 , /?e^/ 2 )) (68) 

with <fi = 2\t. The states \S±) are examples of two-modes cat states of the 
even and odd type and may be characterized by oscillatory distribution of the 
vibrational quanta and by strong nonclassical correlations between the modes. It 
is important to underline the possibility of detecting the entanglement between 
the modes, which is reflected in the interference between the two components of 
|5+). This can be experimentally realized in the same way as in the one mode 
case, that is by measuring the probability P_ (</>) of finding the ion in its internal 
ground state. Indeed P—((f>) is given as (with a and (5 real for simplicity) 

P_(0) = (S+\S+) = \ {l + e -(« 3 +/3 3 )(l-™sW) cog [( a 2 _ ^ gin0 ] | (6Q) 

and is sensitive both to the sum and to the difference of the squares of the 
amplitudes of the coherent states. 



b ) Mesoscopic superpositions of pair coherent states 

Till now we have been dealing with quite direct generalization of the one 
mode cat state ([35]) to two dimensions. There has been however much inter- 
est also in the superpositions of macroscopically distinguishable bimodal states 
such as pair coherent states, SU(2) coherent states, squeezed states and so on. 
These superpositions are usually referred to as Schrodinger cat-like states. Gou, 



Steinbach and Knight [ 155 1 present a scheme for generating mesoscopic super- 



positions of two pair coherent states separated in phase by </>: 

\bq,4>)=Nt{\Z\q)+e i *\-ti;q)) (70) 



In equation (|70| ) is a normalization constant and the pair coherent states 
\£;q) are defined by equation (|59"|). It is easy to verify that the states \£;q, <fi) 
are eigenstates of the operator (a x a y ) 2 with eigenvalue £ 2 . The scheme for gen- 
crating such states is a generalization of the one proposed by the same authors 
for engineering pair coherent states and discussed in the previous paragraph. 
Also in this case the desired target state is obtained as steady state solution of 
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a master equation describing the effect of spontaneous emissions on a system 
represented in terms of the following effective interaction Hamiltonian in the 
interaction picture: 



(a x a y ) - £ 2 



(71) 



For general initial states of motion the oscillatory steady state, detected by 
the extinction of the fluorescence emitted by the ion, is a statistical mixture. 
However, the authors show that, if the initial state of motion of the ion is 
prepared in the Fock state \q, 0) (\q + 1,1)), then the vibrational steady state is 
described by an even (odd) pair coherent state. 



c) Mesoscopic superpositions of SU(2) coherent states 

Let us now consider some recent proposals for generating Schrodinger cat 
like states based on conditional measurements of the internal state of the ion. 
The first scheme we wish to describe is aimed at creating superpositions of 



two SU(2) coherent states separated in phase by n 129, 156|. The system 
considered has been already described in subsection III B 3. We have there 
seen that a single bidimensionally and isotropically confined ion, driven by an 
appropriate configuration of laser beams (see equation(pl|)) and prepared in an 
initial vibrational Fock state along the direction x with an angle ir/4 relative to 
the x axis, exhibits a temporal behaviour strongly sensitive to the parity of the 
initial total number of vibrational quanta. It is shown in ]129[ that, if the ion is 
initially prepared in the vibrational state \n s = N, n y = 0) and in its electronic 
ground state, then the system dynamics, governed by the Hamiltonian (pl[), 
is characterized by the existence of instants of time at which the conditional 
measurement of the electronic ground state projects the cm. motion into a 
superposition of two macroscopically distinguishable SU(2) coherent states. In 
more detail, if the initial total number of vibrational quanta N 3> 1 is even, 
then at t = t e /2 = ttN/AH' the following superposition is generated 

\i> (t e /2)) = ^=(\t= 1,3 = N/2) + (-l)f |t = -l,j = N/2)) (72) 

whereas, if N is odd, the state generated at t — t a ~ t e /2 is 

m<>)) = (|r = i, j = N/2) - *(-1)^|t = —i, j = N/2)) (73) 

Equations (Q) and (|73| ) sh ow th at the properties of the state generated following 
the procedure envisaged in |l29| , strongly depend on the parity of the total initial 
number N of vibrational quanta. In other words this means that, for even N, 
the state (%•)), defined by equation (|72|), has the form of an even (odd) SU(2) 



coherent state [157] , if ^ is even (odd). On the other hand, for odd N, the 



two states \ip(t )), obtained measuring at t = t the internal state of the ion as 
|— ), may be called SU(2) Yurke-Stoler like coherent states, with a difference of 
J (%) in the relative quantum phase, when is even (odd). 

It is in teresting to look at the states ( |72| ) and (|73| ) under another point of 
view Il28|. 
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The angular momentum L z — 2T1J2 can be express in terms of the right and 
left circular quanta operators, n. r = a\a r and hi — aJa/ respectively, as follows 



[158 



L z = h (a\,a r — ajdij (74) 



It is not difficult to verify that the eigenstates \ni = N,n r — 0) and \ni = 0, n r — 
N) of L z , corresponding to its minimum (— KN) and maximum (+HN) eigenval- 
ues respectively, coincide with the SU(2) coherent states |r — i,j — N/2) and 
|t = — i,j = N/2). Thus, the two contributing terms in the quantum super- 
position (|7^) are eigenstates of the orbital angular momentum correspondent 
to its minimum and maximum eigenvalue. In Fig.^ (a) we report the spatial 
distribution of the vibrational state ([73]) generated with the help of our scheme. 
This plot clearly shows interference effects which represent a signature of the 
quantum nature of this superposition. Note that the spatial interference fringes 



are sensitive to vibrational decoherence |15£] and can, thus, be used to study 
the decoherence induced transition of the quantum superpositions ([73]) into the 
correspondent statistical mixture whose spatial distribution is shown in Fig. ^3] 
(b). The vibrational state @, generated in case of even initial N, is on the 
contrary, a superposition of two states describing ionic oscillations along the two 
orthogonal directions x and y respectively. 



e) ID Schrodinger cats of bidimensionally confined ions 
We conclude this section considering an example of the procedure for gen- 
erating ID cat states exploiting laser beams along the x — y plane. Usually in 
these schemes the ion is supposed to be confined in a highly anisotropic two- 
dimensional trap. To generate a given state along the x axis of the trap, it is 
assumed that the ion is tightly bound in the y direction and that the laser a long 



this axis is tuned to the atomic transition. The proposal we wish to review [ 160 1 
deals with the generation of even and odd squeezed states of the form: 

\<P ± )K(\a,$±\-a,0) (75) 

where \a, £) is the squeezed state defined by \a, £) = D(a)S(£,)\0) , with D(a) = 
exp{aa) — a*a) and5(£) = exp^a) 2 — £*a 2 ) displacem ent an d "squeeze" operators 



respectively. Cat like states as the one of equation ( |3.3.3[ ) may be engineered 
as dark state solutions in presence of dissipation due to spontaneous emission. 
Indeed it is not difficult to prove that even and odd squeezed states satisfy an 
eigenvalue equation of the form 

(- gi e-^a 2 x - g 2 e-*4 2 _ 2 g e-^aia x ) p = (p (76) 

In order to ensure that the steady state solution effectively approaches a even 
or odd squeezed state of the form ( 3.3.3j ), how ever , we need to assume that the 



initial state is prepared with a precise parity [ 160 1 . For example starting from 



the Fock state |0) or |1) one obtains even or odd squeezed states respectively. 

4 Conclusive remarks 

Hamiltonian models describing the coupling between pseudo-spin and bosonic 
dynamical variables are ubiquitous in physics. The reason is that when one is 
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faced with a complex problem like the interaction of a material subsystem (atom, 
molecule, crystal, etc.) with the radiation field (electromagnetic or elastic), such 
models catch the essential ingredients of the physical problem and often are ex- 
actly solvable. In this paper we have dealt with a special class of such models, 
namely those describing the coupling of a few-level atom with a pair of bosonic 
modes. We have shown that these models are not oversimplified representations 
of given physical situations only, but that, in fact, most of them have been real- 
ized in the CQED and trapped ion experimental realms. Our review discusses 
a wild collection of results achieved over the last few years in both contexts 
elucidating that their occurrence is chiefly due to the entanglement which gets 
established in the tripartite system. In particular the possibility of inducing 
entanglement dependent quantum correlation effects in the two modes system 
is shown. This very attractive aspect as well as the generation of many nonclas- 
sical bimodal states, clearly witness the dynamical richness of the hamiltonian 
models which we have focused on in this review paper. 
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Figure 1: Energy level diagram of a three- level atom (a) in the A configuration; 
(b) in the S configuration and (c) in the V configuration. 

Figure 2: Average photon number in mode 1 for g\ = g 2 = g and n = 20 

Figure 3: Energy level diagram of a three- level atom in the A configuration. 
The detuning A is assumed large enough with respect to E3 — E\ in order to 
adiabatically eliminate the level E 2 . 

Figure 4: Energy level diagram of a three-level atom interacting with the two 
modes of the field. 

Figure 5: Time evolution of the /j, = 1 average photon number, for n = r 2 = 
s = 0, corresponding to the initial condition (a) n — 20 and (b) n = 21. 

Figure 6: Schematic view of the axial section of a Paul trap 

Figure 7: Raman excitation scheme coupling the electronic states \+)and\—) 
and the vibrational levels through a third upper electronic level \s) 

Figure 8: Vibrational Raman scheme leaving the electronic state unchanged 

Figure 9: Scheme for the measurement of the electronic state of the ion by 
means of quantum jumps techniques 

Figure 10: Time evolution of (Ji(t)) for initial total number of vibrational 
quanta N — 20 (gray line) and N — 21 (black line) 



Figure 1 1 : Probability distributions in the process of generating the pair coher- 
ent state \a(3,q = 2 > pcs from the initial coherent state \a,(3 > with \a\ 2 = 2.5 
and |/3 1 2 = 1.5. The phase difference of the two driving beams is (f> — ir, and 
X /(n Lx - SlLy) = 0.0049 

Figure 12: Spatial distribution of the superposition 

-1= (\m = N,n r = 0) +i\ni = Q,n r = N)), for N = 21, against x' = | 

and y' = l with j3 = y^p 



Figure 13: Spatial distribution of the statistical mixture of the two states |n; = 
N, n r = 0) and \m = 0, n r = N), for N = 21, against x' = ^ and y' = ^ with 
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